Brane inflation and the WMAP data: a Bayesian 
analysis 



Larissa Lorenz 

E-mail: loreiiz@iap.fr 

Institut d'Astrophysique de Paris, UMR 7095-CNRS, Universite Pierre et Marie 
Curie, 98bis boulevard Arago, 75014 Paris, France 

Jerome Martin 

E-mail: jmartin@iap.fr 

Institut d'Astrophysique de Paris, UMR 7095-CNRS, Universite Pierre et Marie 
Curie, 98bis boulevard Arago, 75014 Paris, France 

Christophe Ringeval 

E-mail: ringeval@fyma.ucl.ac.be 

Theoretical and Mathematical Physics Group, Center for Particle Physics and 
Phenomenology, Louvain University, 2 Chcmin du Cyclotron, 1348 Louvain-la-Neuve, 
Belgium 

Abstract. The Wilkinson Microwave Anisotropy Probe (WMAP) constraints on 
string inspired "brane inflation" are investigated. Here, the inflaton field is interpreted 
as the distance between two branes placed in a flux-enriched background geometry 
and has a Dirac-Born-Infcld (DBI) kinetic term. Our method relies on an exact 
numerical integration of the inflationary power spectra coupled to a Markov-Chain 
Monte-Carlo exploration of the parameter space. This analysis is valid for any 
perturbative value of the string coupling constant and of the string length, and includes 
a phenomenological modelling of the reheating era to describe the post-inflationary 
evolution. It is found that the data favour a scenario where inflation stops by violation 
of the slow-roll conditions well before brane annihilation, rather than by tachyonic 
instability. Concerning the background geometry, it is established that \ogv > —10 at 
95% confidence level (CL), where v is the dimensionlcss ratio of the five-dimensional 
sub-manifold at the base of the six-dimensional warped conifold geometry to the 
volume of the unit five-sphere. The reheating energy scale remains poorly constrained, 
Trch > 20GeV at 95% CL, for an extreme equation of state (wrch ^ —1/3) only. 
Assuming the string length is known, the favoured values of the string coupling and 
of the Ramond-Ramond total background charge appear to be correlated. Finally, the 
stochastic regime (without and with volume effects) is studied using a perturbative 
treatment of the Langevin equation. The validity of such an approximate scheme is 
discussed and shown to be too limited for a full characterisation of the quantum effects. 



PACS numbers: 98.80.Cq, 98.70.Vc 
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1. Introduction 

Inflation is presently considered our best description of tlie early Universe, since it 
solves the puzzles of the Standard Hot Big Bang Model [U [21 El IH O [6] and, in addition, 
provides a very convincing mechanism for structure formation leading to an almost 
scale- invariant power spectrum [71 El [9l [10] . The small deviations from scale-invariance 
are connected to the micro-physics of inflation and are, therefore, of utmost importance 
when investigating the physical conditions in the very early Universe [HI [121 llH]- In 
practice, inflation is usually driven by a scalar field (or possibly several fields) and 
for this type of matter, the effective pressure can become negative, if the scalar field's 
potential is sufficiently flat. This condition is necessary to produce a phase of accelerated 
expansion (i.e. inflation) in General Relativity. In this context, one of the main open 
issues is the physical nature of the inflaton field. Since it is clear that this question 
should be addressed in the framework of extensions of the Standard Model of particle 
physics. String Theory seems to be a promising place to look for a physically well- 
motivated candidate field. In recent years, this issue has given rise to several studies, 
for reviews see references [HI [151 llHl [111 [HI [19]. In particular, string inspired brane 
world scenarios have turned out to be especially fruitful. In this context, the inflaton 
field is usually interpreted as the distance between two branes moving relative to each 
other in the extra dimensions [20l [2T1 [22] . The corresponding potential generically 
contains a Coulomb term of the form ~ C — where C is a constant and the inflaton 
field [221 [231 [21 [251 [261 [271 [281 [29] . The end of inflation can then either occur by violation 
of the slow-roll conditions or by a mechanism of tachyonic instability, as in hybrid 
inflation. In either case, when the distance between the branes becomes of the order 
of the string length, the simple single field description breaks down and a more refined 
(fully stringy) treatment is necessary. The main goal of the present article is to carry out 
a detailed comparison of the inflationary predictions derived from one particular class 
of brane inflation scenarios, namely the Kachru-Kallosh-Linde-Maldacena-McAllister- 
Trivedi (KKLMMT) models, with the Cosmic Microwave Background (CMB) data 
presently available. 

Analysing this class of potentials is further motivated by two rather different 
reasons. Firstly, although it is possible (for instance by using the slow-roll 
approximation) to span the space of (single field) inflationary models without relying 
on a specific form of the inflaton potential [301 [211 [321 [33] , it is also interesting to not 
use any such approximation except, of course, the linear theory for the cosmo logical 
perturbations [34l [35l [36]. Let us also mention that with the incoming flow of more 
accurate data from the Planck satellite [371 13H] or from the future arc-minute resolution 
CMB experiments [39l [^0] [H] , the errors linked with the current approximation schemes 
may become a problem. The no-approximation approach often requires numerical 
computations. As a consequence, one has to pay the price of choosing a particular 
form for the inflaton potential, which implies some loss of generality. A possible way 
out and a systematic strategy to span the space of inflationary models is to choose 
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potentials that are representative examples of a whole class of scenarios. For instance, a 
scenario where inflation takes place in the regime of large (in Planck mass units) vacuum 
expectation value (vev) for the field, and where the potential goes to infinity as the vev 
of the field increases towards infinity, is illustrated by chaotic and/or hybrid inflation. 
The difference between these two models consists in the way inflation stops. On the 
contrary, a model where inflation occurs for small values of the inflaton vev is well 
represented by models in which the potential has a symmetry breaking shape. These 
classes of models were studied in great detail and compared to the WMAP third-year 
data in reference [31]. From this general perspective, the KKLMMT models with their 
above-mentioned Coulomb term generate inflation when the vev of the field is large, but 
with a potential bounded at infinity. 

Another motivation, different in spirit from the previous one, lies in model building 
issues. The Coulomb type potential naturally arises in string inspired scenarios of 
inflation. Moreover, this type of potential has been known for a long time, even without 
relying on stringy inspiration, since it can occur in models of hybrid inflation from 
dynamical supersymmetry breaking |32]. Recently, in the context of String Theory, 
corrections to the Coulomb potential associated with the problem of moduli stabilisation 
have been explored [TTJ [2ll [251 [SSI [27], [29] . Although interesting from the theoretical 
point of view, they are not considered in the following. Including these effects in the 
analysis would greatly increase the number of free parameters relevant for the description 
of the inflationary part of the model and their degeneracy by opening the parameter 
space. As a result, given the current data, no constraint could be obtained (as it is 
the case, for instance, for running mass inflation [M]). We have therefore chosen to 
perform an exhaustive Bayesian analysis of the CMB data by considering only the 
generic Coulomb contribution to the potential in the context of string inspired brane 
inflation. 

Compared to the existing literature |l3l [lH [15], our study uses a complete 
Markov-Chain Monte-Carlo (MCMC) analysis based on an exact numerical integration 
of this type of potential. In this respect, we do not consider particular values of the 
string length y/a' and the string coupling gg but, on the contrary, develop a parameter 
scanning strategy which allows to leave these parameters fre^. However, we will see in 
the following that, already with only the Coulomb part of the potential, degeneracies 
in the parameter space do not permit the WMAP data to fully constrain these stringy 
quantities. Nevertheless, the throat geometry and some combinations of the KKLMMT 
model parameters are constrained. Finally, we present new approximate solutions for 
the field evolution in the DEI regime and refine the analysis of inflationary quantum 
effects. In particular, it is shown that the approximate solutions usually considered are 
not sufficient to fully characterise this regime. 

This article is organised as follows. In section [2], we briefly recall some basic facts 
about the KKLMMT scenario of brane inflation, deriving the governing equations from 

I Some conditions on these parameters must be nevertheless fulfilled as, for instance, < 1 in order 
to be in the perturbative regime. 
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the effective four-dimensional action based on type IIB String Theory. In particular, 
we discuss the intrinsic string features of the model, which among others manifest 
themselves through an unusual kinetic behaviour as well as a restricted domain of 
validity for the model. In section [3l we briefly consider the impact of quantum 
fluctuations that may affect the classical trajectory in certain regions. Section H] is 
devoted to the slow-roll phase of the scenario, which can be treated analogously to small 
field models. However, in the case of the KKLMMT scenario, the slow-roll discussion 
has to be complemented by a study of the stringy aspects: this is presented in section |5l 
After deriving the restrictions on the parameter space from model-intrinsic arguments 
in section [HI we present our results from the WMAP third-year data in sections [7] and [HI 
Conclusions are drawn in the final section [91 Finally, in Appendix A , we give details on 
the stochastic regime. 



2. Brane inflation 



In the following, we discuss the KKLMMT model of brane inflation, first proposed in 
reference |22] as a realisation of inflation within String Theory. Our goal here is to 
consider all aspects (string related or not) at play during the inflationary phase. We 
begin with a intuitive description of how the KKLMMT model arises in String Theory. 



2.1. Qualitative description 

The model considers a D3 and an anti-D3 brane in a ten-dimensional supergravity 
background whose world-volume (time-axis included) is aligned with the coordinate 
axes, /i varying from to 3. The six extra dimensions with A = 4, . . . , 9, are 
compactified such that the "radial" distance separating brane and anti-brane is y"^ = r, 
and their distance in the other coordinates vanishes. While inflationary models involving 
D3 (or higher dimensional) branes have been considered in the literature before, the 
KKLMMT model assumes that the six-dimensional section forms a so-called Klebanov- 
Strassler (KS) throat [171 [HI [IHl [501 EI], i-e. a warped deformed conifold with 
background fluxes for certain background fields. 

Intuitively, this background can be thought of as populated by A/" heavy D3 branes 
with Ramond-Ramond (RR) charge. The anti-D3 brane is embedded at a fixed position 
To (the bottom of the KS throat), representing an additional source of RR field strength 
that can be described as a small perturbation. The D3 brane is inserted into this 
perturbed background at ri ^ rg (far from the bottom of the throat) and is considered 
as a probe: it does not affect the geometry itself, but experiences the forces due to 
gravity and RR interaction, brought about by the exchange of light closed string modes 
between the branes. The distance r = ri — ro between the test D3 and the fixed anti-D3 
brane corresponds (up to normalisation) to the inflaton field whose potential V{(j)) 
can be calculated from the forces experienced by the test D3 brane in the limit r ^ is, 
is = y/a' being the string length scale. Inflation takes place while "rolls down" along a 
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flat interaction potential in the direction of decreasing 0. When the branes become too 
close, in the sense that their proper distance approaches the string scale, a tachyon, the 
lightest open string mode stretching from one brane to the other, appears, and the long 
distance potential V{(j)) is no longer valid. When the D3 brane reaches the bottom of 
the throat ro, the two branes annihilate in a complex process whose details are beyond 
the scope of the present paper. We will, however, introduce a phenomenological model 
which assumes that the brane annihilation triggers a reheating era. 



2.2. The inflaton effective action 



In this section, we use an effective field theory representation of the KKLMMT model, 
and our starting point will be the effective action of the inflaton field 0. Considering 
that the inflaton cj) is itself an open string mode, it is related to the distance r between 
the brane and the anti-brane by 

0=V^r, (1) 

and its effective four- dimensional action reads 



S 



where the quantity k is defined by 
Svr 




,2 ' 



and 



1 



(2) 



(3) 



(4) 



is the D3 brane tension, being the string coupling constant. The quantity ^(0) is 
the (perturbed) warp factor of the ten-dimensional metric whose form can be calculated 
from the full Einstein equations [5^. To simplify the notation, one can introduce the 
effective brane tension T{(j)) 



(5) 



allowing us to re-write equation ([2]) as [52j 
S 

1 

~ ~2k 



-g d X. 



(6) 



It is clear that, for a slowly varying scalar field, the terms in (90)^ are small and the 
square root in 5'dbi can be Taylor expanded in the field derivatives. The resulting action 
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at leading order identifies with the one of a canonically normahsed scalar field evolving 
in a self-interaction potential. 

2.3. The kinetic term 

On a flat Friedmann-Lemaitre-Robertson- Walker (FLRW) brane, the metric reads 

ds^ = a(r/)2 (-dr/2 + Sijdx'dx^) , (7) 

where r) and are the conformal time and spatial coordinates. In terms of the cosmic 
time dt = adr], the metric tensor reduces to diag [— 1, a^(t), a^(t), a^(t)]. For such a 
homogeneous metric, can only depend on time and keeping the first non-trivial term 
in the square root expansion reproduces a standard kinetic term, namely 



a%t)d^x, (8) 



where a dot denotes a derivative with respect to the cosmic time. However, this remains 
true as long as higher order terms in the expansion can be neglected, i.e. for 

0'<T(0). (9) 

Let us denote by 0dbi the field value from which onwards the field evolution would 
be dominated by the string corrections to the standard kinetic term. Notice that 
equation Qj is also the condition required to have a positive square root argument 
in equation ([6]). In this sense, ■\/T(0) can be understood as the field maximum speed in 
the warped throat, and one can, by relativistic analogy, define a Lorentz factor 7 such 
that [53] 

7(0,0) = ^ ^ (10) 

VI - <pyT (0) 

When 7 ~ 0(1), the field (p follows the dynamics of a canonically normalised scalar 
field which can slow-roll (and hence, produce inflation) on a sufficiently fiat potential. 
However, there also exists an "ultra- relativistic" regime, 7^1, where the string- 
intrinsic DBI form of the kinetic term is important. In this regime, inflation may be 
possible even if the scalar field potential is not flat. 

2.4. The inflaton potential 

Having established the inflaton dynamics, we now identify the potential V{(j)). There is 
only one term in the expansion of equation ([6]) which does not contain field derivatives 
and hence we conclude that 

V(0) = 2T(0) = ^ . (11) 
h{(j)) 

So far we have not specified the form of the warp factor h{(f)). In the case of a KS 
background perturbed by one anti-D3 brane at the bottom of the throat, we have 
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Figure 1. Left panel: The dotted blue line represents the potential given by 
equation (|13p . The solid green line shows the approximate potential for 3> /i, see 
equation (|14p . The conventional slow-roll phase occurs while the field is rolling on the 
extremely flat region to the right of the plot. Right panel: Same as the left panel, but 
in logarithmic unit for the potential. 



leading to the potential 



Vi^) = ,^. , ,,4 • (13) 
1 + 



In the case where ^ /i, the first order expansion in reads 

4" 

V{^) ^ 1 - ' " ' 



(14) 



which is the most commonly used form of the KKLMMT Coulomb-type potential. The 
full expression of the potential as well as its approximate expression are represented in 
figure [H 

2.5. Parameter sets 

From the String Theory point of view, a complete description of the KKLMMT model 
is given by the parameter set 

{g,,a\v,M,IC). (15) 

where we remind that and y/a' are the string coupling and length, respectively. 
The quantity v is the (dimensionless) volume ratio of the five-dimensional sub-manifold 
forming the basis of the six-dimensional conifold geometry to the volume of the five- 
sphere 

Vol(X.) 

For instance, the base of the KS throat is the Einstein space Ti i, hence in this case one 
has V = 16/27 [51]. The quantity A4 and /C are two positive integers associated with 
the separately quantised background fiuxes. They satisfy 

J\f = )CM, (17) 
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where A/" is a positive integer representing the total background RR charge. 

Another set of parameters equivalent to f|T5l) is given in terms of the extra 
dimensions' geometry. The two natural length scales are the location of the bottom 
of the throat tq, which corresponds to the field value 

00 = V^ro, (18) 

and the throat edge ruv The position r = ruv where the throat is smoothly glued into 
the rest of the extra dimensional bulk. Provided the depth of the throat is comparable 
to its width, the edge field value can be approximated by [23] 

0UV = Vr^r^v, (19) 

with 

r^, = Arrg^a"-^ . (20) 

Consequently, one can deduce the maximum warp factor between the edge and the 
bottom of the throat [55] 



/ -^^P hn^ ' (21) 

where the anti-D3's small perturbation to the KS warp factor has been neglected. The 
set of geometrical parameters is therefore given by 

ig,,a',(Po,<Pvv,^f). (22) 

The potential expression in equation ( fT3il contains the parameters M and /i, which 
are defined from the string parameters as 

Note that v can be expressed through M and /i. 

Since both M and /i are proportional to the ratio (p^/M , one has to provide another 
parameter to maintain equivalence with the set of equations (|T5l) and (|22|) . for instance 
M. In that case, equation ( l23i) can be inverted for 



'2t^ \ V 

The model is thus equally well described by the "cosmological" parameter set 

{g,,a\M,^l,M). (26) 

2. 6. Domain of validity 

We move on to discuss the domain of validity of the effective field theory description. 
It turns out that in order for equation ([2]) to be justified, one has to impose some 
constraints on the model parameters and the accessible field values 0. 
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2.6.1. Throat within the overall six- dimensional volume. In equation ([3]), the four- 
dimensional gravitational coupling constant k is expressed as a function of the four- 
dimensional Planck mass, which, in a ten-dimensional geometry, depends on the volume 
of the compactified extra dimensions, 

T^tot 

ml, = 87T-^, (27) 

where Kiq is related to the string parameters by 

/tio = ^{27cyg^,a'\ (28) 

and Vq°^ is the total volume of the six compactified extra dimensions [56l [57]. Now 
remembering that the total six-dimensional volume consists in the volume V^*'^''"'^* of 
the KS throat, where inflation takes place, plus the bulk volume V^*'""^, into which this 
throat is glued at r = r^y, one has 

ytot ^ ybulk ^ ^^throat^ ^29) 

implying V^°yV^^'°^^ > 1. Although we ignore the precise form of the six-dimensional 
bulk, the throat volume can be calculated from [21] 

^^throat ^ 2n%^^a"rl,, (30) 
with Tuv given by equation ( |T9l) . Then, re- writing equation f l271) as 

T^throat T/tot 

= 8^ ; ^tLoat > (31) 

gives the total to throat volume ratio 



V / V 

^ = «'<y4vr3^s^. (32) 

Since this ratio must exceed one, we obtain a condition on the possible combinations of 
J\f and V, depending on the choice of gs and a': 

_^3/2^-i/2 ^ 27r'/'a'ml,gl/'. (33) 

Note that this condition heavily depends on the choice of gs and a', the value of which 
can vary significantly according to the specific String Theory considered. It is also 
convenient to recast equation ( l33l) in terms of 0uv- From equations (fT9l) and (!32i) . one 
has 

A T ^throat 

-*Sv = ^^. (34) 
and the volume ratio limit fl33|) implies 

Let us also remark that, for the supergravity effective description to be valid, the volume 
of the extra dimensions, and therefore of the throat, should be larger than a'^, even when 
the volume of the five- dimensional conifold basis is very small. However, one can show 
that this condition does not play a role in the following analysis. 
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2.6.2. Inflation within a single throat. We will restrict our attention to the case where 
the D3 brane begins its journey towards the anti-D3 brane inside the KS throat, i.e. for 



where 0uv is the field value on the throat edge [see equation (fTOll ]. In cosmological 
terms, this means that all the e-folds required to solve the Standard Big Bang Model 
problems occur inside one and the same throat. We will discuss in section |6] in which 
sense equation fl36l) allows us to impose restrictions on the parameter values. 

2. 6. 3. End of brane motion by an instability-like mechanism. The derivation of f|T3l) 
assumes that brane interaction is due to the exchange of the lightest closed string modes 
only. This is only approximately true: Firstly, one could include the contribution of 
heavier modes (whose propagation through the bulk is Yukawa suppressed) as well 
as consider more general effects such as the Kahler potential necessary to stabilise 
the moduli. In this case, deriving the full brane-anti-brane interaction potential is 
highly non-trivial. A phenomenological approach is to assume that these effects can be 
summarised in an additional potential term [e.g. of the form V{(j)) = m?(fp'] which can 
be added to equation (fT3ll . This route has been explored in the literature, see e.g. ^5] . 
Very recently, the theoretical underpinning of form and origin of these terms has also 
received much attention [241 EHl EHl [27l [28] . For our cosmological purposes, these effects 
have not been considered, following the argumentation outlined earlier based on the 
accuracy of the present data. Secondly, in addition to the above-mentioned corrections 
to the potential, the effective field model intrinsically has a restricted domain of validity. 
When the two branes come too close in terms of their proper distance s (with respect 
to the string scale 4), an open string can stretch between them and a tachyon, being 
the lightest excitation of this open string, appears. Let us determine the value 0strg for 
which the proper distance between the branes Send — 4, where 4 = V^- It is solution 
of the following equation [15] 



where the origin sq can be chosen at the brane annihilation location tq, i.e. Sq = 0, 
with To the radial coordinate of the bottom of the throat. Using the warp factor of 
the unperturbed background, h{r) = r^y/r'^, one gets Send = '"uv ln(rstrg/'ro)- Since 
Send = V^, we have rstrg = roe^/^^v^ which after normalisation gives 



where use has been made of equations (fT9l) and ( l25l) . Notice that 0strg exphcitly depends 
on the background fiux J\f. As a consequence, another equivalent formulation of the 
parameter set ( l26l) is therefore given by 



< 0, 



(36) 




(37) 




(38) 



(fl's,a',M,/i,0strg)- 



(39) 
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From equation (1351) . one immediately gets 

0strg > /W, (40) 

since the background flux number A/" > 1. Although the String Theory details of the 
tachyon appearance and brane annihilation are out of the scope of this work, we will 
assume that these processes trigger a reheating era which precedes the cosmological 
radiation-dominated era (see section IHTTl) . 

2. 7. Equations of motion 

The equations of motion for the infiaton in a FLRW metric on the brane can be obtained 
from the total action ([6]). From the stress-tensor associated with the scalar field 0, the 
energy density and pressure read 

p = (7 - 1)T(0) + y(0), P = ^T(0) - 1^(0), (41) 

7 

from which one obtains the Friedmann-Lemaitre equations 

H' = ^[T (0) {^-1) + V (0)] , H = ^T{<P) - 7) • (42) 

When the velocity of the field is small, the Lorentz factor 7 defined in f fTOj) is close to 
one and its square root can be expanded. Likewise, varying the action with respect to 
the infiaton field leads to the following "Klein-Gordon like" equation of motion 

d^0 , /3 , d\nH\ d0 , ^3 ^\ r{ct>) , V\4>)-T\ct>) _^ .... 

dAT^ + ^72 ^ dN ) AN^ V7' 7 2/72 ^ 
where, in the present context, a prime denotes a derivative with respect to 0. Note that 
A^ here is the number of e-folds (not to be confused with the background flux M) , which 
will be used as a convenient measure of time. In the limit 7^1, equations (1421) and 
( 143|) reduce to the standard Friedmann-Lemaitre and Klein-Gordon equations. It will 
be important to establish in which regime the DBI corrections to the kinetic term (and 
therefore, deviation from the standard slow-roll dynamics) may be important. 

2.8. Field evolution overview 

Brane inflation according to the KKLMMT scenario incorporates both features of usual 
slow-roll inflation as well as intrinsic String Theory effects. Let us sketch the successive 
phases of the field evolution as rolls down the potential V{(f)) towards its small values. 

In figure [2] (left panel), one can distinguish a region (for very large field values, dark 
green shaded) in which the quantum character of the infiaton field cannot be neglected: 
Quantum fluctuations are of the same order as the classical ones where the potential 
is extremely flat, i.e. for exceeding a certain 0fiuct • This is notably of interest if the 
brane starts its journey above this limit. The peculiar properties of this regime will be 



discussed in section |3] and Appendix A 



For the moment, let us assume that the brane motion starts at 0in < 0fluct- In 
this region, the potential (fT3ll is still very flat since ^ p, allowing for most of the 
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Figure 2. Left panel: Sketch of the expected dynamical regimes according to the vev 
of the inflaton field for the potential ([T3)) . The field starts out in the flat (no hatched) 
region of the potential and rolls towards smaller field values. Eventually, it will reach 
the light green hatched region where the slope of the potential becomes noticeable and 
where the DBI effects can no longer be neglected. This regime may or may not be 
reached according to the value of ipstrg for which the derivation of equation (jl3[) breaks 
down. On the far right (dark green hatched region), the potential is extremely flat and 
quantum fluctuations are expected to dominate over the field classical evolution. Right 
panel: The slow-roll parameters ei (green line) and £2 (blue line) for the KKLMMT 
potential. The solid curves have been obtained for /i/mp, = 0.1 whereas the dotted 
ones for n/mpi = 0.01. Decreasing /i/mpi increases the differences between and 
(j)^^, the field values for which ei = 1 and £2 = 1, respectively. Notice that we always 
have 4>e2 > 0ci ■ 



inflationary expansion to take place. One may therefore expect the usual slow-roll 
approximation to be valid and it will be used in section H] to derive the shape of the 
induced primordial power spectra. 

Rolling further to the left, the field eventually enters a region where the 
slope of the potential becomes noticeable (light green shaded region). Therefore, the 
conventional slow-roll approximations are certainly no longer sufficient to describe this 
evolution, and, since the field velocity should increase, one may expect the DBI dynamics 
to be the driving force for < 0dbi- In addition, 0strg is also located around this potential 
domain. The precise order of events in this region will be discussed in section [51 



3. Stochastic inflation 



In this section, we will denote the classically predicted field value by 0ci, and consider 
the effect of the quantum fiuctuations on the classical trajectory [58l EHJ EOl [61]. The 
field values for which quantum fiuctuations are comparable to the classical ones can be 
obtained by first estimating the size of the classical fiuctuations in the slow-roll regime. 
Denoting V^^j the potential derivative with respect to the field, one gets 

A0ci^-T^^^At, (44) 

3-f/(0cl) 
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which comes from the Klein-Gordon equation by neglecting 0. For At it is convenient 
to take the typical time scale at work during expansion, namely one Hubble time, 
At = 1/H. Concerning the quantum fluctuations, one has for a test field in de-Sitter 
space-time. 



A0qu 

The field value 



27C 



(45) 



Wet 5 for which classical and quantum fluctuations are of equal 
amplitude, verifies A0ci(0fluct) = ^(pqa{4>auct) ■ From equations (IHI) and (145!) it follows 
that 

1/10 

(46) 



?fluct 



/i 



3 

8^ 



M 



Using equation fl23l) . we can express (j^6l) in terms of v and the overall scale of inflation 
M/rupi. According to their value, 0fluct can generically be much larger than the minimal 
initial field value needed to produce A^t — 10^ e-folds of inflation. To ensure the 
consistency of the semi-classical approach, we require 

0in < 0fluct , (47) 

in the following. Quantum effects can be nevertheless described using the stochastic 
approach [60j [61] which is presented in detail in Appendix A Since the initial value of 
the field is also limited by the requirement of having the D3 brane inside the throat, the 
effective semi-classical model remains valid only if < min(0f)uc;t; 0uv)- 



4. The slow-roll regime 

In this section we analyse the KKLMMT potential according to its slow-roll 
characteristics. When denoted as in (fT3|l . the analogy between V^(0) and small field 
models (investigated in detail in reference [34j) is evident. Notice, however, that for 
small field inflation, is increasing while inflation is under way, whereas KKLMMT 
inflation proceeds from large to small field values (see figured]). 



4-1- Slow-roll parameters 

Among the several definitions of the slow-roll parameters, we here use the Hubble-flow 
functions e„ defined by [62], 163] l64] 

dln|e„| , , 

'"+^^^iV ' ^^^^ 
where is the number of e-folds since some initial time r][^. The above hierarchy starts 
from eo = Hi^/H. Slow-roll inflation proceeds as long as |e„| <^ 1, for all n > 0, while 
inflation takes place if the scale factor is accelerating, i.e. ci < 1. 



Brane inflation and the WMAP data 



14 



From the potential fll3p . the first two Hubble flow parameters in their slow-roll 
formulation read 

10 

mp, \ \Ji 



1 

ei = - 

TT 



£2 = - 

TT 



/i 



1 + 



2 ' 



5 + 



/i 



2 ' 



(49) 



(50) 



and they are represented in flgure [2] (right panel). They depend on the ratio 0//i and the 
vev energy scale /x/mpi [or M/nipi, by virtue of equation (123|) ] only. As above-mentioned, 
inflation ends when ei = 1, while the slow-roll approximation itself breaks down when 
€2 = 1 (and/or ei = 1). Let us stress again that the fleld theory description of the brane 
motion fails at = 0strg- 



4-.2. Classical field trajectory 



The next step is to obtain the classical fleld trajectory; this is possible while the slow-roll 
approximation is valid but even in this case, the trajectory is found only implicitly. The 
total number of e-folds N{(j)) reads [65] 

6 i/i\2 i/j.\6 T/j\2" 



N{(j)) = 2n 



'ml [e (, j + 2 (, /i J 6 UJ 2 UJ J ■ ^^^^ 

For the generic situation where ^ 1, the two quadratic terms are sub-dominant. 
Neglecting them yields the following explicit solution 

1 1/6 



0(AO 



3 

TT 



mp 



. (52) 

/i y n \ jjL J 

As seen in flgure |2], given that the fleld evolution starts on the very flat part of the 
potential (but outside the domain of quantum fluctuations), there will surely be an 
initial phase of evolution during which the slow-roll approximation, and hence the 
trajectory fl52l) . is valid. However, moving towards smaller fleld values, this trajectory 
will eventually become inaccurate and one has to determine at which fleld value the 
slow-roll approximation breaks down and how this fleld value compares to the end of 
inflation and/or to 0strg- 



4-3. When does the slow-roll regime break down? 

The slow-roll approximation breaks down when one of the slow-roll parameters becomes 
of order one, which need not coincide with the end of accelerated expansion (ei = 1). 
If slow- roll breaks down flrst, inflation can still proceed (typically during a very small 
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number of e-folds), but the evolution of the field is no longer described by the slow- roll 
trajectory (|52|) . Let us define the field values and by 

ei(0,J = l, e2(0ej = l. (53) 

From equations (l49l) and (!50!) . one obtains the algebraic equations 

and 



7r \ /i 



/i \ / \ 5 / mpi \ / /i 



1 + 2^ + 



1 / \ 2 / \ 10 



TT \ /i 

0. (55) 

n \ fi ' ■ 

These equations cannot be solved explicitly (except numerically). Some approximate 
analytical solutions are derived in the following. 

4-3.1. Case fi > ntp^ with </>ei,e2 < Iii this case, keeping only the two dominant terms 
in equation flMl) gives 



NO 1/ \2 / \10 

\ 1 / mpi \ / /i 



^£1/ ^ V / \V^ei 

whose solution is 

mp, 1 



0, (56) 



^ • (57) 

This expression is consistent with our assumptions: /i > mpi indeed implies < /i. 
Still in the same limit, the three remaining terms in equation (155!) are 



xo r / \2/ \6 1/ \2/ xlO 

jjL \ of mpi \ / \ 1 / '^pi \ / 



-£2/ \ y VC/),^/ 71" V / \re2 

This equation has an exact solution with two acceptable roots 

^£2 \ ^Pl ' 2 



/i / , fi \ n 




. (58) 



(59) 



For /X > mpi, slow-roll inflation can proceed till the field reaches max(0ej,0e2) < A^- 
By comparing equation (!57|) and (!59|) . we further see that > 0^^, i.e. the slow-roll 
approximation breaks down before the end of inflation. 

However, as already shown in section [5], the KKLMMT model itself is no longer 
well-defined for < 0strg due to the stringy origin of the potential f|T3l) . The value 0strg 
given in equation (l38l) is always greater than /x. Hence, in the case where fi > rripi, 
inflation will definitely come to an end at 0strg > niax(0fj, 0^2)- -^^ a result, the entire 
field evolution occurs in the slow-roll regime and, as it will be shown in the following, 
we do not need to worry about DBI effects which remain negligible in that case. 
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4-3.2. Case jj, < nipi with 0e^,e2 > In this limit, the second and third term of 
equation (IMl) are small compared to one, which leads to 



Similarly, in equation ( l55l) . the second and third term are small compared to one, while 
the last term is small compared to the fourth, so by keeping only the two dominant 
terms we obtain 



Comparing equation (!6T!) and (!60|) shows that one still has 0^2 > 4>ei, and the slow- 
roll approximation breaks down before inflation stops. This is confirmed in figure [2j 
However, this time, the field does not necessarily reach (pstrg first and non-slow-roll 
inflation may continue between and 0ei- Moreover, from (p,,-^ to 0strg, the field 
encounters a transitory regime interpolating between the conventional slow-roll and 
DBI dynamics. As discussed above, all the slow-roll formulae are strictly speaking no 
longer valid after e2 = 1. As a result, the region < 0^2 is only accessible numerically 
and we will show in this way that the number of e-folds occurring between 0^2 and 0^^ 
(or 0strg) is actually small. The key lesson to be learned from 0^2 > 4'ti is that whenever 
we want to rely on analytic results of the slow-roll calculus, we have to confine the field 
to values > 0e2 for consistency and keep in mind that the results will be applicable 
only if the number of e-folds occurring after 0^2 remains negligible. 

To proceed further, let us summarise the possible field evolution given that fi < rupi: 
The first option is that, also in this case, the field reaches the value 0strg first, 
i.e. 0sti.g > 0^2; this is possible depending on the values of J\f and v for the chosen 
background geometry [see equation (l38l) ]. In that case, the entire field evolution 
corresponds to slow- roll inflation and is described by the analytical trajectory (l52l) . 

If 0strg < 0e2) then the situation is more complex and depends on the exact order of 
0ei, 0DBI and 0strg- For instance, if (p^i is the next field value to be reached after 0^2, the 
accelerated expansion would come to a halt here for a field with standard dynamics, 
rendering the number of e-folds between 0^2 and 0^^ extremely small. However, the 
subsequent evolution for < 0^^ could bring the field in the regime of DBI dominance 
and inflation may re-start. At the same time, 0strg could be located between 0^^ and 
0DBI and, as a consequence, 0dbi would never be reached. All other combinations are of 
course a priori possible. Which one is realised in practice depends on the value of the 
free parameters characterising the model. 

5. String-intrinsic aspects 

We now move on to investigating those aspects of the KKLMMT model due to its stringy 
origin. In particular, two questions remain to be properly addressed: the unusual DBI 




(60) 




(61) 
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dynamics, and the appearance of a tachyon once the proper brane distance approaches 
the string scale. 

5.1. When does inflation end? 

The potential (fT3l) results from String Theory under the assumption that only the 
lightest closed string modes (gravitons and the RR particles) contribute to the brane 
interaction. When the proper distance s between the branes becomes comparable to 
the string length scale 4? the exchange of heavier closed modes become relevant, and 
an open string can stretch from one brane to the other. The appearance of a tachyon 
triggers the brane annihilation process, or reheating from a cosmological point of view. 

We already calculated the field value 0strg at which the tachyon appears in 
equation fl38l) . The crucial question is, does the field first reach the value </)strg or 
0^2? If < 0strg; inflation ends at brane annihilation for (p^nd = 0strg and the field 
spends its entire "lifetime" in the usual slow-roll regime. Such is the course of events 
when fi > mpi. For /i < mpi, one has to compare (p^^ given in equation (161]) . with (pstrg 
given by equation (l38l) . Their ratio reads 



i^strg 



exp 



(62) 



This expression involves the energy scale M/rrip, of the potential because (p^^ and 0strg 
both depend on M/mpi, but not in the same way. However, from an observational point 
of view, this scale is involved in the amplitude of the cosmological perturbations and 
certainly fixed by the CMB normalisation. Therefore, we may go further and access the 
above ratio using the current WMAP measurement of the CMB quadrupole. 

Let us call 0* the field value at the time when the wavelength of the cosmological 
perturbations of physical interest today left the Hubble radius during inflation. The 
slow-roll field trajectory allows us to express (p^ in terms of A^,,, the number of e-folds 
between the time of Hubble exit and the end of inflation at (pend- Using the approximated 
classical trajectory ( l52i) . this leads to 

t^ ) ~ \ t^ ) ^\t^) 
To proceed further, we have to specify when exactly inflation ends. But, at the same 
time, this is also the question we try to address. Changing the mechanism which stops 
inflation will change the scale M/iripi and, hence, will affect the ratio ( l62i) . 

It is convenient to determine the frontier 0^2 = 0strg in the parameter space. Setting 
for convenience 0end = 0^2 equation (163|) gives 

1/6 

_ • (64) 



mpA / 5 



TT \ /i / \ 3 

Note that, had we used 0end = 0ei instead, the result would be fl64|) up to the replacement 
A"^ + 5/3 ^ A"^. Therefore, ending inflation at 0^2 instead of 0^^ only causes a small 
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shift in N^. This is expected since, as soon as 62 > 1, slow-roll is violated and the field 
starts to evolve rapidly: even if the precise evolution can only be probed numerically (see 
section [521) , the number of e- folds spent in this regime is generically small. Inserting the 
result (!64|) into the expressions (149!) and (!50!) for the slow-roll parameters, one arrives 
at 



ei ~ - 

TT 



rrir- 



£2 



3N, 



1 + 



3N, 



rup 
-1 



-5/3 



(65) 



(66) 



As can be seen in these equations, the value of ei is generically much smaller than €2- 
Moreover, the quantity €2 does not depend on the scale n/rupi, at first order. The CMB 
quadrupole normalisation gives the additional relation 

rms— PS 



where, for 0=,, ^ fi, one may approximate K 
ei, we find the relations 

M 



(67) 

M^. Using the above expression fl65l) for 
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The numerical value of the WMAP quadrupole is 

Qrms-PS 



T 



6 X 10" 



(68) 
(69) 

(70) 



Inserting this expression into equation ( l62i) . the ratio ^ej/^strg becomes a function of Af 
and f only, at fixed string coupling Qs- Explicitly, one obtains 



-'stre 



exp 



47r5fs— 

V 



-1/4 



where the constant C reads 



C = 10'/%6N, + 10 f I 45^^^:^ 



-1/8 



(71) 



(72) 



In the plane {v, A/"), the condition 



domains 0strg < 



and 



•'strg 



> 



= 0strg is a curve separating the plane into the two 
,. Notice that the shape and position of the frontier 
depend on the value of gs] this is illustrated in the left panel of figure [31 

It is convenient to re-scale the parameters Qs, M and v to minimise the dependence 
on of the contour (j)^^ = (pstrg- For this purpose, we define the new variables 



X 



V 



V = 



(73) 



Usually, the analysis is performed with some specific values of gg- In this paper, our 
strategy is different and more general. With the help of the above rescaling, our results 
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Figure 3. Left panel: The contour (f)^^ — 0strg in the plane (In t), In A/"), obtained from 
equations ((62)) and ([7T|) using the normalisation given by the CMB quadrupole with 
N^: = 50 [see equation ([68]) ]. The dotted line corresponds to gs = 0.1, the dashed line 
to gs = 10^'' and the dotted-dashed one to gs — IQ^^ . In each case, the area enclosed 
by the contour is the part of the parameter space where the slow-roll conditions are 
violated before brane annihilation and (/)strg does not play an important role for the end 
of inflation. It is clear from this plot that the contour sensitively depends on the value 
of gs- However, such a dependence can be absorbed by an appropriate rescaling of the 
parameters as shown in the right panel. The same contour (jj^^ = 0strg is represented 
in the plane (In a:, Inw), these parameters being defined in equation ([73]). It is universal 
for all values of the string coupling gs- 



will be valid for any values of the coupling constant. In terms of these new parameters, 
the condition (p^trg = 4'e2 from equation (ITTl) now reads, in logarithmic units 

Inw = 81nC - 21nx - 8a;"^/^ (74) 

This contour is represented in figure [3] (right panel). All dependence of this line on 
Qs has been absorbed in the rescaling, so one can now state universally that inside 
the contour, slow-roll breaks down before the model-intrinsic instability is reached 
(0strg < 0e2)) whereas outside, slow- roll inflation proceeds all the way until brane 
annihilation (0strg > <Pe2)- fact, this result is more general and we demonstrate in the 
following that the rescaling permits to absorb the Qs dependence in all the equations 
expressing a physically relevant condition for the model. 

Note that in the above discussion, we have used 0end = to derive the 
normalisation fl68l) . In doing so, we have in fact ignored the other "stringy" characteristic 
of the KKLMMT model, namely the DBI kinetic term. While brane evolution definitely 
ends at ^strg, it is possible that, due to the unusual dynamics, a phase of "DBI inflation" 
occurs even after violation of the slow-roll conditions. If a considerable number of e-folds 
could be produced in the DBI regime, the value of used in equation (|68ll would no 
longer be correct. In the following, we discuss for which field value (pY^Bi those effects 
will be important and how much expansion may occur outside the slow-roll regime. 
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5.2. When is the DBI regime important? 

An order of magnitude of the field value 0dbi for which the DBI regime is relevant can 
be obtained from the potential (fTTi) by use of the condition ([9]): 

^ ^(0db:) . (75) 

For a standard kinetic term (7 = 1), this value would precisely coincide with ei = 1. 
This suggest that the DBI regime appears when standard inflation ends. Of course, given 
the increasingly non-canonical dynamics as the field is approaching 0dbu the standard 
slow-roll formula is no longer valid and should be replaced by its generalised form 

dlnH 1 fd<py 2 fdlnHV 

where use has been made of equation fllOj) and of the exact expression Therefore, 
in the ultra- relativistic limit 7^1, equation fl76l) suggests that the expansion of 
the universe may accelerate even with a steep potential. If such a second phase of 
inflation occurs for more than typically 60 e-folds, then the slow-roll phase becomes 
observationally irrelevant and the DBI regime crucial. We address this issue by analysing 
the DBI regime using both analytical and numerical methods. 

The above expression for ei can be further simplified by using equation ffTOl) to 
express the field derivatives in terms of 7, H and T solely. The Friedmann-Lemaitre 
equation (142|) can then be used to remove any explicit dependency in the Hubble 
parameter in favour of 7: 

H' = K . (77) 

1 + 7" 

As a result, the first slow-roll parameter reads, in terms of 7, V and T 
3 1-7'' 



Equations (177|) and (178!) are exact and clearly enhance the effect of 7. In the limit 
7^1, one recovers the usual slow- roll formulae obtained with a standard kinetic term, 
whereas an approximate analytic solution can be derived in the ultra-relativistic limit 
7 3> 1. Under the assumption V/T = 0(1), at leading order in 7"^, one gets 

61 = ^ + 0(7-^), (79) 

showing that deep in the DBI regime the universe is not inflating but expands as in a 
matter dominated era. In our case V/T = 2 which ensures the validity of the above 
expansion. Notice, however, that this relationship does no longer hold if one considers an 
additional term in the potential, e.g. m^0^. In that case the V/T term in equation (1781) 
can no longer be neglected and inflation may indeed proceed in the DBI regime B6]. 
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Although the universe is not inflating for 7 ^ 1, we still have to check that the 
number of e-folds the field spends in the DBI "matter dominated" era remains small. 
From equations ( 1761) and ( 1791) . the Hubble trajectory at leading order in 7"^ reads 



H{N) = Hi exp 



2 ^ 



Ni) 



(80) 



where Hi and A^^i respectively denote the Hubble parameter and the e-fold at which 
the limit 7^1 becomes relevant. They could, for instance, be defined as the ones 
corresponding to the end of the slow-rolling phase provided that the transitory regime 
occurring between the conventional slow-roll and the region of DBI dominance is short. 
The field evolution at leading order is given in terms of 7 by equation (!76i) and reads 

2 



The system of equations is closed by using expression fH2l) for the Hubble parameter 



H' 



1) 



which is exact since, in our case, one has V/T 
trajectory (neglecting one compared to 7 ^ 1 

d(f) ^ 
dN ~ 



Hi 



exp 



12) 



2. As a consequence, one gets the field 



^3) 



This expression can be implicitly integrated in terms of the Gauss hyper-geometric 
function through use of equations (ITT]) and (IT^ . yielding 
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where x = fi/(f). A more illuminating form can be obtained in both limits x ^ 1 and 
X <C 1. In the limit x -C 1, that is to say ^ /i, using the Taylor series defining the 
hyper-geometric function, equation ( |84l) becomes 
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which gives for the field 



exp 
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(86) 



Similarly, in the limit x ^ 1 or <^ /i, using the linear transformation formulae 
associated with the Gauss hyper-geometric function [67], equation fl84|) yields 

3V2H1 jj 
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A^i + - In 
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Figure 4. Evolution of the field </>, the Hubble- flow functions ei, £2 and the DBI 
parameter 7 in the last e-foldings of an extreme model living close to the throat edge 
(/i = mpi/V327r). This model has been chosen to emphasise the DBI effects: the 
solid lines correspond to an exact numerical integration of the full action ^ whereas 
the dashed lines are obtained by using the slow-roll approximations (|50p and 

([5T|) . The DBI regime smoothly connects the slow- roll evolution (7 ~ 1) to the ultra- 
relativistic matter-like expansion (7 ^ 1). Note that in the model at hand, brane 
annihilation occurs at 0strg > M preventing any observable effects coming from the 
DBI evolution. 



and, therefore, inverting the previous relation, one obtains 

-1 

exp - (A^ - A^i) - 1 ' 



3V2H1 



(88) 



2 MM 

In fact, the hmit x ^ 1 (0 <^ /i) cannot be reached in the model under scrutiny since 
the tachyonic instability occurs for (p ~ (/)strg > A*- As a result, the ultra-relativistic DBI 
regime could only occur for (p > fi. As can be seen in equation fISB]) . the corresponding 
field evolution is exponentially fast implying that the number of e-folds during which 
(p > 0strg and 7 ^ 1 is negligible, provided 4>i/fi and the denominator are not too 
large. This is indeed the case for the following reasons. Firstly, as already mentioned, 
equation fl78|) shows that 7 can deviate from unity only if ei is also of order unity, i.e. the 
field should be in the non-flat region of the potential and therefore 0i > fi. Concerning 
the denominator, as can be seen from equation (1771) . Hi/M"^ < l/'^pi and this term is at 
most of order fi/rripi. The constraint (l35l) coming from the size of the throat, together 
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with the fact that brane annihilation occurs at 0strg > 1^, require that 



m 



PI 



(89) 




ensuring that Hifi/M'^ cannot be large. 

To end this section, we have numerically checked that the previous analysis was 
still qualitatively valid during the intermediate regime in which neither the slow-roll nor 
the ultra-relativistic approximations can be used. Figure H] shows the last e-foldings of 
evolution for an extreme model, close to the limit (|89l) . for which the effects of the DBI 
regime can be seen. The differences between the slow-roll approximation (under the 
standard kinetic term hypothesis) and the exact DBI integration appear only during 
less than an e-fold and close to the region where ei = 1. Remembering that the effective 
field description of the brane motion physically ends at 0strg, we conclude that for all 
practical purposes, the DBI regime has no cosmological observable effects for the model 
based on the potential (|T3l) . The situation, however, can be entirely different in models 
that include additional terms (e.g. m^0^) in the potential, as shown in reference |45j . 

6. Theoretical restrictions on the parameter space 

We are ready to conclude our analytical investigation of the KKLMMT inflationary 
model by deriving restrictions on its free parameters. For our numerical calculations and 
for comparison with the WMAP3 data, we use the observable parameter set (M, /i. A/"), 
where M and /i are related by equation (12^ . There are essentially three theoretical 
consistency relations to be satisfied. 

6.1. Constraints from the size of the throat 

For the model to be consistent, the warped throat in which the D3 and anti-D3 branes are 
located should be smaller than the total size of the compactified sub-manifold spanned 
by the extra dimensions. Earlier, we derived the restriction (!33l) from this condition. 
It turns out that this condition can be conveniently re-written in terms of our rescaled 
parameters (x, v) defined in equation ( |73l) as 



This is a straight line in the plane (Inx, Inv) whose offset depends on a' (see figure [5]). 
A significant fraction of the parameter space is cut-out by this requirement. Notice 
that, as announced above, the rescaling ( |73l) has removed any dependence of this bound 
on Qs and it is universal in the sense that it does not involve the precise inflationary 
trajectory, nor the mechanism that ends inflation. 

6.2. Constraints from starting inflation inside the throat 

For all the e-folds to occur inside one throat, we enforce < 0uv- Since decreases 
during inflation, this restriction applies to the initial field value 0in. On the other hand. 



Inv < In (a'mpjVr) Inx. 



(90) 



Brane inflation and the WMAP data 



24 




n(47Tg^/v) ln(4Tig3^/v) 



Figure 5. Allowed regions (ticks in) in the rescaled parameter plane (lnx,ln-y) for 
a'mpj — 1000 (left panel) and a'rn^^ = 10 (right panel) with the fiducial values N^, ~ 50 
and = 60. The black dotted curve is the contour (p^^ — (pstrg of figure [3l For all 
points located inside this contour, the slow-roll conditions are violated before brane 
annihilation. The green dashed line (with ticks down) represents the volume ratio 
constraint ([90)) . whose slope is universal but whose offset depends on a'mp,. Regions 
above this line are therefore excluded. The solid green curve (with ticks up) represents 
the condition that all of the brane evolution occurs within one throat, and has been 
obtained through a numerical integration. All points below that curve would not 
satisfy this condition. Its shape can be piecewise analysed. In the region 0strg < 4>il2 
(inside the black dotted contour), this is a straight line given by equation ([9^ . The 
slope of this line is universal, but the offset again depends on a'rn^^. Outside the 
dotted contour, namely for 0strg > 0£2i and in the limit 0strg ^ </'e2J this boundary 
is described by equation (jlOip . Combined, these pieces lead to the solid green curve 
with ticks up; to make the shape of both pieces visible individually, they have been 
extended outside their respective domains of validity (dotted-dashed black curve) . 




there is a lower bound on 0in such that at least iVx e-folds are produced. This bound in 
turn depends on the mechanism ending the inflationary expansion. 

Let us start with the case in which the slow-roll approximation is violated before 
brane annihilation. Making use of the classical trajectory (!52|) . the initial field value 
necessary to produce A^t e-folds before can be derived from equation (!6TI) and reads 

3 + iVxJ I ■ (91) 

For the KKLMMT case, imposing the CMB normalisation ( l69l) calculated earlier and 
using equation (fTOl) . one can recast the bound 0in,e2 < 0uv in terms of the rescaled 
variables 

Int; > -Inx -41nX), (92) 
where the constant T) reads 

V = {Tia'miy'^' (6iV, + 10)^/^^ (6iV, + 10)-^/^ (^45^^^^ . (93) 
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This is also a straight hne in the plane (x, v) and further constrains from below the 
parameter space (see figure [5]). As before, all the gg dependence has been removed and 
only the a' one remains in the expression of the offset. This offset also depends on A^^, 
and A^^ but only weakly because these two quantities appear in a logarithm. In the 
following, the fiducial values A^* = 50 and A^^ = 60 will be used 

Notice that by using the normalisation fl69|) derived from 0end = (pt2^ ^^e previous 
result is valid only inside the contour fl74l) . Therefore, we still have to treat the case 
where inflation ends prematurely at 0strg > <Pe2- Since the entire field evolution occurs in 
the slow-roll regime, the trajectory (!52|) remains valid all the time. The minimal initial 
field value 0in,strg leading to A't e-folds of inflation is therefore obtained by setting 
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where 0strg is given by equation ( l38l) . In order to use the CMB normalisation, the new 
value of 0=1, has to be derived from equation 



/i 



^strg 



3 

TT 



/i 



-2 



A^. 



where, now, 0end 

1/6 



•'strg- 



One obtains 
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The slow-roll parameter associated with this field value reads 
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and can be plugged into equation (1671) to give the CMB normalisation condition 
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Unlike the case 0end = 0e2 5 ^^is equation for M/rUpi cannot be made explicit but one 
can use equation flM]) to replace the last term in the square brackets and obtain 
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An explicit analytical solution can be derived in the limit 0strg ^ 4>e2- Due to the 
presence of the A^,,, term, this approximation may be violated only for large A^* values. 
But if we are not in this extreme situation, in the limit 0strg ^ 0e2) gets 
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Figure 6. Plots similar to those in figure [5l for a'm^^ — 1000 on the left and for 
a'nipi = 6.4 (as chosen in [22]) on the right. The quantities and are fixed to 
their fiducial values iV, = 50 and A^^, — 60. The additional symbols correspond to 
various models selected according to the values oi v, Af and gg. All purple symbols 
correspond to the KS throat value oiv — 16/27, and respectively have A/" = 1, 10, 100, 
10'^, 10^, 10^ and 10^ from left to right. The diamond-shaped points are for gs = 0.1, 
the triangles for = 10"'^, and the squares for g^ — 10^^. The original KKLMMT 
model [22] is indicated by a red cross in the right panel; it corresponds to w = 1, 
TV = 160 and gs = 0.1 at a'ml^ = 6.4. 



+ 2\nv + ^\nx-3\n{7ia'ml,) - In (^45^^^^^ . (101) 

Again, there is no trace of Qs left after rescaling, but instead of a straight hne as was the 
case for equations fl90l) and fl92|) . this is an imphcit contour in the (Inx, Inv) plane. It 
is represented in figure O Let us recall that equation (1 10 II) is only valid far enough from 
the contour 0strg = such a way that 0strg ^ 0e2- ^^e intermediate regime where 

inflation does end by instability at (pstrg, but the two terms in the square brackets of 
equation (|98|) are of comparable value, the correct normalisation of M/rripi to the CMB 
quadrupole is only accessible numerically and has been also represented in figure [5] for 
convenience. It is remarkable that when one normalises with either fl68|) (obtained from 
0end = 4>t2) with ([99]) (assumiug (pstrg ^ (p€2)y the gs dependence can be absorbed into 
the unique rescaling given in equation (175]) . 

Finally, in figure [6l various concrete models for different values of v, J\f and Qs 
are compared with the consistency conditions derived above. In particular, one notices 
that the volume ratio constraint is quite restrictive and that many possible models are 
already ruled out by this condition. 

6. 3. Constraints from stochastic inflation 

If the field (p starts out at a value comparable to </)fluct given in equation (H6|) . then 
stochastic effects will be important and the classical field trajectory does no longer 
suffice to describe the field evolution. As a minimum requirement, in the case where 
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brane annihilation occurs after (p^^j has to impose 

4>m,e2 < 0fluct, (102) 



where 0in,e2 is given by equation flQlj) . Comparing equation fl9T|) with equation fH6|) . 
and using the correct normahsation ( l68l) . we notice that the dependence on A/" and w 
disappears. As a result, the restriction coming from (11021) concerns the maximal total 
number of e-folds A^t achievable, 

12 \ -3/5 



N.<-1+ (322")-'/' (6A^, + 10) (^45^?^^ 



which gives max(A^T) — 10^ for A^,, ^ 50. Similarly, if inflation ends at (j)stTg > 4>e2J 
imposes 

<^in,strg < </>fluct, (103) 

to avoid quantum fluctuation dominance. In this case, using equations (H6l) and 
together with the normalisation given in equation fl99l) . gives 

In l + 6iV^ (45 '^™-''" ^ '"i-'"ii-"3exp(-5i-i/i 



<5ta2Jl„(45%^). (104) 

As usual, the rescaling allows us to get rid of any explicit dependence in gs. We do 
not pursue this issue further here. For a detailed study of the presence of quantum 



fluctuation effects in the throat, see [Appendix A 



7. Implications of the WMAP slow-roll bounds 

In this section, the theoretically predicted values of the KKLMMT slow-roll parameters 
are confronted with the third year WMAP data (WMAP3). The WMAP3 bounds on 
the Hubble-flow parameters ei and €2 have been discussed in reference [Ml [331 EH] under 
various prior choices and hypotheses. For the sake of robustness, in the following we 
use the constraints derived in reference [31] on the first order Hubble-flow parameters 
obtained by marginalising over the second order slow-roll parameter €3 and under an 
uniform prior choice for log(ei) in [—5,0] and for €2 in [—0.2,0.2]. The resulting one 
and two-sigma contour intervals of the two-dimensional marginalised posteriors are 
represented in figure [71 

Let us first consider the situation where violation of the slow-roll conditions occurs 
before brane annihilation. In this case, if one chooses for instance A^* = 50, the CMB 
normalisation from equations (l68l) and (l69ll yields 

— ~ 5.6 X 10"%"^/^ — ~ 1.4 X 10"^t;"^/®. (105) 
nipi rupi 

Hence these two parameters become functions of v only. The resulting numerical values 

for ei and 62 are given by equations f[65l) and f[66l) and read 

ei -4.8 X 10"i^t;"^/^ €2 ~ 0.03 . (106) 
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This has two very important consequences: it is clear that except for very small values of 
the volume ration v < 10^^*^, ei remains negligible compared to €2- In particular, since 
the spectral indices of the scalar and tensor primordial power spectra, at first order in 
the Hubble-flow parameters, read 

ng - 1 = - 2ei - €2, = -2ei, (107) 

we have unobservably small tensor modes and a scalar spectral index 

Us ~ 0.97. (108) 

In figure [71 we have studied the slow-roll predictions in more detail. The slow-roll 
parameters have been calculated exactly in the sense that equation ( |55l) has been solved 
numerically, that is to say we have obtained the exact value of 0^2- The same has been 
done for 0^,. Then, the corresponding values of the Hubble- flow parameters have been 
derived with the help of equations ( l65l) and ( l66l) for different N^, such that A^* G [40, 60]. 
For each value of A^,,,, this gives a point in the space (61,62) and for the range under 
consideration, a segment line spanning twenty e- folds (see figure [7]). The "left" end of 
this segment line corresponds to the largest value of A^^,, namely A^^, = 60 and the "right" 
end corresponds to A^* =40. The left panel assumes v = 1, while the right one was 
plotted for v = 10~®. The dotted-dashed blue contours are the WMAP3 constraints 
on the slow-roll parameters and corresponds to the 68% and 95% confidence intervals 
of the two-dimensional marginalised posteriors. The dotted black lines are the lines of 
constant ng, the red one tracing the scale invariant case ng = 1. The figure confirms the 
previous analysis. The gravitational waves level is generically very low and the spectral 
index is slightly red in full agreement with the WMAP3 data. Moreover, changing the 
size of the extra dimension changes ei without affecting 62. However, for reasonable 
values of w, ei remains so small that there is no hope to detect primordial gravitational 
waves in the future. 

Let us now turn to the situation where the parameters of the model are such that 
inflation ends by instability at 0strg- In this case, the WMAP normalisation leads to 
the scale M/rripi given by equation ( !99l) . Then, straightforward calculations lead to the 
following expressions for the Hubble-fiow parameters 
7;V3 / n2 \ -V3 / , N -20/3 

^ I .^Vrms-PS \ I Vstrg \ 



57;l/3 / n2 \ -1/3 / , N -8/3 



.^^^45— J ^-.j . (110) 

In the limit 0strg//^ ^ 1, one has ci ^ £2 ^ 1 except maybe for extreme values of 
V. Therefore, rig is now pushed towards one, while becomes even smaller for all 
values of A^,,,. As a consequence, one expects this scenario to be disfavoured by the data. 
Indeed, as can be noticed in figure [3, a scale invariant power spectrum is not likely if 
the gravitational waves level is low (this is no longer true if ei > 10~^). 

To conclude this section, let us recap what can be deduced from the slow-roll 
analysis. Firstly, the situation where inflation ends by violation of the slow-roll 
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Figure 7. Left panel: Predictions for the KKLMMT model with w = 1 in the case 
where one considers that inflation stops at 62 = 1 . Right panel: Same as left panel but 
wfth V ^ 10"^. 



condition is favoured by the data as compared to the case where inflation stops by brane 
annihilation. Secondly, one expects the level of gravitational waves to be extremely low. 
The stringy interpretation of this result is linked to the so-called Lyth bound [691 El] 
which relates the tensor to scalar ratio to the total variation of the infiaton field value. In 
our case, this gives A0 = a/Ts (r^^ — tq) — V^r^^ = Therefore, given equation 
f l35|) . the volume bound on the throat limits the gravitational wave contribution, as 
recovered here. Thirdly, the constraint ci < 3 x 10~^ implies the limit 

log (47r2w) > -16 . (Ill) 

This is clearly not a very stringent constraint. In the next section, we will see that, 
combined with the constraint on the size of the throat, the data can lead to a better 
limit. 

8. WMAP constraints in the exact numerical approach 

In this section, the cosmological consequences of the KKLMMT model are investigated 
using a numerical integration of the brane motion up to its linear perturbations to extract 
the exact scalar and tensor primordial power spectra. The tachyonic instability occurring 
at the bottom of the throat triggering a reheating era is considered through a simple 
phenomenological model. The predicted CMB power spectra can then be computed by 
integrating the seeded cosmological perturbations through the radiation and matter era 
and we have used a modified version of the CAMB code for this purpose [70]. This method 
allows a Markov-Chain Monte-Carlo (MCMC) analysis of the WMAP third year data 
involving the usual cosmological parameters together with the KKLMMT parameters. 
A modified version of the COSMOMC code [71] has been used to derive the probability 
density distributions of the cosmological, KKLMMT and reheating parameter values 
given the data. 
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8.1. Method and hypotheses 

The numerical method used has been introduced in references [34l [35] and consists in the 
mode by mode integration of both the background and the perturbed quantities till the 
end of inflation (see also Refs. [72l[73l[71l[l5l[76l[77l[78]). As shown in section [E2l the 
DBI regime is not relevant for the model under scrutiny and we have chosen to integrate 
the background equations stemming from the action ([6]) in the standard kinetic-term 
limit. The resulting equations of motion are obtained from equations ( H3l) . ( 1761) and 
([77D in the 7 = 1 limit and for the KKLMMT potential (fT3D. 

In the longitudinal gauge, we consider the scalar and tensor linear perturbations to 
the flat FLRW metric (CD, 

ds^ = {-(1 + 2(^)dr]'^ + [(1 - 2^)6ij + hij] dx'dx^} , (112) 

where $ and \Ef are the Bardeen potentials [79] and hij is a transverse and traceless 
tensor. From the perturbed Einstein and Klein-Gordon equations, the dynamics of the 
scalar and tensor modes, in Fourier space, reduces to the equations of motion of two 
uncoupled parametric oscillators [TJ Eni EH ES] 

^ + .;3%(fc,7^KT = 0, (113) 
dr] 

where the frequencies are respectively 



u;Uk,ri) = k'-^-^^, ujl{k,r^) = e--, (114) 



a" 



the prime denoting differentiation with respect to conformal time. The scalar and tensor 
mode functions are gauge invariant and read [7] 

/is = aV2K I 50+ — — 1 , = ahij , (115) 

where stands for the linear perturbations, in the longitudinal gauge, of the infiaton 
field [see equation ([1])]. The initial conditions for the perturbed quantities /Xg X are 
set in the decoupling limit k ^ 7i where the deep sub-Hubble modes behave as free 
quantum fields. Assuming their initial state to be the usual Bunch-Davies vacuum, one 
gets [83] 

Q-ik{v-Vi) 

lim /is,T = 7= — , (116) 

with rji some (arbitrary) initial conformal time. These initial conditions uniquely 
determine the solutions of equation (11131) which can be numerically integrated along 
the lines described in reference [35]. From such a numerical integration, the primordial 
scalar and tensor power spectra can be evaluated at the end of inflation for all the 
relevant observable wavenumbers 



= k^P^{k) 



V, ^ ep,{k) = ^ ^ . (117) 
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Concerning the background quantities, the attractor mechanism during inflation 
ensures that the initial fleld value has no direct observable effects since it determines 
only the total number of e-folds A^t- Its value can therefore be, a priori, arbitrarily 
chosen provided Nrj, is big enough to solve the problems of the standard hot Big-Bang 
phase and allows the observable perturbations to be sub- Hubble initially. As thoroughly 
discussed in sections [5] and O the brane motion ends by the appearance of a tachyon 
from its coalescence with the anti-brane at the bottom of the throat, and this determines 
the fleld value at which the above classical evolution cannot longer be used, namely 
~ 0strg- The initial field value 0in can therefore be chosen in such a way that the brane 
motion generates A^t e-folds of inflation before reaching (pstrg- However, one still has to 
impose that the brane motion starts inside the throat for the model to be consistent, 

i.e. 0in < 0UV- 

According to the previous discussion, it is natural to assume that the reheating 
era starts when the brane alights on the anti-brane. Although, strictly speaking, brane 
annihilation occurs at the bottom of the throat ro, the String Theory details of the 
evolution for (p < 0strg are out of the scope of our simple approach [HI [851 186] . 
As a result, we have adopted the phenomenological reheating model discussed in 
references [Ml [35] . It assumes that the reheating can only influence the observed 
perturbations through its effects on the cosmological redshift -Zstrg, which will be 
identifled with the beginning of the reheating era. Such an assumption is motivated 
by the fact that adiabatic super-Hubble perturbation modes should not be signiflcantly 
modifled during reheating, at least in the absence of entropy modes (here, we deal with a 
single fleld model). The influence of 2;strg can be seen by rewriting equation flll3p in terms 
of the number of e-folds as a time variable. The wavenumbers to be considered during 
inflation always appear in the ratio k/Ti. which should be determined from the observable 
wavenumbers k/aQ measured today, typically three decades around k^: = 0.05 Mpc^^. 
At a given e-fold A^ during inflation, 

k k Nt-N 

:t7 = -(1 + W^- (118) 

Assuming instantaneous transitions between the successive expansion eras, the redshift 
2;strg is given by 

1 , / o \ , '^strg 1 

where preh is the total energy density at the end of the reheating era, and f^rad, Hq are 
the density parameter of radiation and the Hubble parameter today. The flrst two terms 
in equation (I119p are clearly reheating dependent and in absence of a microscopic model 
we can deflne a phenomenological reheating parameter R^ad such that 

InR.ad^-ilnf^'j+ln^. (120) 

4 VPstrg/ O-reh 

This parameter has a simple physical interpretation: it encodes the global deviation the 
dynamics of the reheating era may have with respect to a pure radiation-like era (for 



ln(l + Zstrg) = 7 In (/tVrch) - In - 7 In (\^3Q^Ho) , (119) 

4 OrpVi Z \ / 
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which Rrad vanishes). Notice that, from the numerical integration, the energy density 
Pstrg when (f) = {/)strg is known and only depends on 0strg and the potential parameters. 

Under the previous assumptions, the KKLMMT inflation era requires the knowledge 
of five primordial parameters: the KKLMMT potential parameters, namely M and /i, 
the reheating parameter Rrad and the field values 0uv and 0strg encoding the observable 
properties of the throat and the branes. The resulting CMB anisotropies are uniquely 
determined once the "low energy" cosmological model is fixed. We are considering 
the ACDM fiat universe model which adds four cosmological parameters: the number 
density of baryons f2b, of cold dark matter Qdm, the reduced Hubble parameter today 
h, and the redshift of reionisation Zre- 

8.2. MCMC parameters and priors 

Due to parameter degeneracies with respect to the CMB temperature and polarisation 
angular power spectra, some parameter combinations are more appropriate for an 
efficient MCMC exploration. Concerning the ACDM parameters, instead of directly 
sampling models according to the values of the set (fib, ^dm, ^re, -f^o), it is more 
convenient to use the equivalent set (fib^^, ^dm^^, t, 6'), where r is the optical depth 
and 9 measures the ratio of the sound horizon to the angular diameter distance [TT] . 

Similarly, it is more convenient to use an optimal derived set for the primordial 
parameters. The amplitude of the scalar primordial power spectrum at a given 
observable wavenumber = Vq{K) is a well measured quantity [see equation (|671) ]. 
It is therefore more convenient to directly sample the models according to the values of 
rather than the potential energy scale M/mpi. This can be done by integrating the 
perturbations with the artificial value M/m^i = 1 and then perform a rescaling of M/mpi 
from unity to its physical value that would be associated with the wanted P*. As shown 
in reference |3l], under the rescaling M/rupi — >■ sM/rupi, the power spectra become 
V{k) sV (s^/^/c). As a result, the value of s required is the ratio V^/V^"^\ where 
V^^"^^ is the amplitude of the scalar power spectrum obtained with M/rupi = 1 and 
evaluated at = k^s^^^"^. In fact, to circumvent such a rescaling on the wavenumbers, 
it is more convenient to introduce the rescaled reheating parameter R defined by 

In R = In Rrad + ^ In (/tVstrg) • (121) 

The quantity R represents an effective energy, in Planck units, at the time of brane 
merging, which is exactly equal to A/^Pstrg a radiation-like reheating. The advantage 
of R with respect to Rrad is that, at fixed R, one has fco = K [31]. Finally, as shown in 
section [5], one may expect some observable effects coming from the situations in which 
there is violation of the slow-roll conditions before brane annihilation, i.e. according to 
the value of 0strg/0e2- Since is determined by /i/mpi only [see equation (|6T1) ]. we have 
chosen to perform the MCMC exploration on the parameter 0strg//^ instead of 0strg- 

Finally, the nine MCMC parameters are: f2b/i^, fidm^^, t, on the cosmological 
side, /i/mpi, 0uv, 0strg//^ for the primordial parameters and R for the reheating era. 
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We still have to specify their prior probability distributions according to the theoretical 
constraints and our best knowledge on their value. For the base cosmological parameters 
fib^^, ^dm^^, T and 9, wide top hat uniform distributions have been chosen centred over 
their preferred value from the previous analysis of the CMB data [871 E!- Fo^^ the 
primordial parameters, V^. fixes the amplitude of the cosmological perturbations and 
we have chosen a uniform prior on the logarithm compatible with the amplitude of the 
CMB fluctuations: 

2.7 < In (10^°P*) < 4.0. (122) 

The parameter /i/mpi is related to the underlying String Theory model by equation (!23ll 
and we will assume that the order of magnitude of the volume ratio v is not known. As 
a result, it is natural to also use a uniform prior probability distribution on log (/x/mpi) 
to ensure its conformal invariance. The same considerations hold for the two other 
primordial parameters, 0uv and 0strg, whose dependency with respect to the fundamental 
string theory parameters is given in equations ([T9l) and (!38|) . Our ignorance on the 
values of A/", Qs and a' motivates the use of uninformative uniform logarithmic priors. 
However, the analysis of section [6] imposes various bounds. The maximal size of the 
throat in equation fl35l) and the requirement Af > 1 give the upper limit 

log(v^0uv) <log2. (123) 

Similarly, the definition (138|) of 0strg ensures that 0strg > implying the lower limit 

log(0,trg//i) > . (124) 

Since 0strg < 4>uv, these conditions impose 

log(v^/i)<log2, (125) 

which will be our upper prior limit for the parameter fi/rripi. However, we should 
ensure that there are enough e-folds of inflation to solve the flatness problem and to set 
the sub-Hubble initial conditions (11161) for the observable perturbations. As shown in 
reference [HE], only for some extreme reheating models the values of may exceed 10^. 
In order to include all the models, we have implemented a "hard prior" which, during 
the MCMC exploration, rejects any model that does not support at least 110 e-folds of 
inflation inside the throat. This is implemented in the following way. Once 0strg is known 
[from log(0strg//^) and \og{y/K^)], one can determine 0in such that there are 110 e-folds 
of inflation in between: the model is accepted if 0in < 0uv and rejected otherwise. The 
prior limits on In R are determined by requiring that the end of reheating occurs before 
nucleosynthesis, characterised by pnuc, and that the instantaneous equation of state 
parameter Wreh satisfies the strong and dominant energy conditions — 1/3 < Wj-ch < 1- 
Under these assumptions, one gets [M] 



^ In (^Vnuc) < In R < In (^Vnuc) (^^Vstrg) • (126) 

This equation clearly involves pstrg and the upper prior limit on In R depends on the 
other primordial parameters. We have therefore coded another "hard prior" checked 
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during the MCMC exploration to dismiss any model violating this limit. The energy 
density pnuc at nucleosynthesis time has been quite extremely set around the MeV scale: 
In R > —46. Finally, we have to set a lower limit on the \og{y/kii) prior. According 
to the values of w, fi/rupi may be extremely small. Our choice has been motivated by 
numerical convenience and we have chosen log^y/Kfx) > —3. Indeed, it turns out that, 
for smaller values, one runs into tricky numerical difficulties. The dependency of the 
results with respect to this choice will be carefully discussed in the following. Notice that 
the lower limit of the log(A/K0uv) prior is implicitly set by the others since 0uv > (pm 
and we ensure that there are at least 110 e-folds of inflation in between and 0strg- 

8.3. Data used 

As already mentioned, the CMB measurements used are the WMAP third year 
data [m [871 IMl EOl- The degeneracies between some cosmological parameters have 
been reduced by adding the Hubble Space Telescope (HST) constraint Hq = 72 ± 
8km/s/Mpc [91] and a uniform top hat prior on the age of the universe between 10 Gyrs 
and 20 Gyrs [7T]. The likelihood estimator is provided by the WMAP tean|§| and we 
have used the current 2.2.2 version. Use has been made of the eigenvalue compression 
option at low multipoles to decrease the computing time. The convergence of the Markov 
chain simulations has been monitored from the Gelman and Rubin i?-test implemented 
in COSMOMC [92j . The MCMC exploration has been stopped once the ratio of the variance 
of the means to the mean of the variances between the different chains was less than 
R — 1 < 0.2%. The total number of samples obtained from this convergence criterion 
and drawing the posterior distributions is around 1.5 millions. In the following, we 
present and interpret the posterior marginalised probability distributions for the nine 
model parameters obtained from the previous priors and data. 

8.4. Cosmological parameters 

The marginalised posteriors and mean likelihoods of the cosmological parameters and 
power spectra amplitude are the same than those found in previous analysis of the same 
data but that were using other primordial power spectra, such as the first and second 
order slow-roll spectra in reference [31] or the phenomenological power law models in 
reference [87]. This is expected since these parameters are rather well constrained and 
supported by the slow-roll analysis of section [7] showing that the KKLMMT model can 
indeed produce the preferred value of the power spectra amplitude and spectral indices. 
The best fit model is found with a — 3538.1 for nine parameters (or 3517 degrees of 
freedom after eigenvalue compression). This may be compared with the WMAP power 
law best fit model having — 3540.8 with three parameters less. Although it shows 
that the KKLMMT model can provide a good fit to the data, it is not currently favoured 
because of its intrinsic numbers of parameters. 

§ http://lambda.gsfc.nasa.gov 
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Figure 8. Marginalised posterior probability distributions (solid lines) and mean 
likelihoods (dotted lines) for the sampled primordial parameters of the ACDM- 
KKLMMT model. 



8. 5. Primordial parameters 

The marginalised posteriors and mean likelihoods for the sampled primordial parameters 
log(A/K/i), log(A/K(/)uv), log(0strg/yu) and In R are represented in figure [HI The fact that 
the mean likelihood (dotted curves) remains uniform for log(v^/i) and log(v^0uv) 
shows that, on the prior range explored, these parameters do not help to improve 
the fit to the data. On the other hand, their probability distributions are not flat 
and show that large values of yu/mpi and small values of 0uv are strongly disfavoured: 
log(y^yu) < —1.1 at 95% confidence level. As discussed in reference [93], discrepancies 
between the marginalised probability and the mean likelihood occur when "volume 
effects" are induced by strong correlations between the parameters. The full likelihood 
in the multi-dimensional parameter space can be uniform along peculiar regions whose 
volume may be parameter dependent. The shape of the yu/mpi and 0uv probabilities 
come from the amount of fine tuning required to have a successful model of infiation 
when yu/mpi is close to the throat edge. In that case, there is indeed not too much 
parameter space available for the brane motion since 0uv has to saturate its maximal 
value to allow both (p^trg > and </)strg < (pin < (f>vv- As a result of these correlations. 
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the marginalised probabilities for fi/rupi and 0uv penalise the values for which strong 
fine-tunings are required to satisfy the theoretical priors. Let us mention again that 
these priors are imposed by the consistency of the model. 

On the other hand, the marginalised posteriors and mean likelihoods associated 
with R and 0strg//^ show that these parameters are directly constrained by the data. 
Clearly, 0strg ^ A* is disfavoured and at 95% confidence level one has log(0strg/yu) < 1-4. 
This result can be understood from section [3 Large values of 0strg/ A* correspond to 
models in which brane annihilation occurs well inside the slow-roll regime, i.e. 0strg > • 
As a result, the observable perturbation modes are generated during the brane motion in 
the fiat part of the potential (see figure Ej) and both the slow-roll parameters are small: 
the spectral index is close to rig — 1 while the amplitude of the tensor modes is negligible 
with respect to the scalar modes. As already mentioned in section [7l this situation 
is disfavoured by the WMAP data. The upper bound on 0strg/At is in fact slightly 
dependant on the lower limit of the fj, prior due to the presence of some correlations 
with the other parameters. As discussed in the following, a //-prior independent upper 
limit can be more conveniently given for the rescaled parameter 0strg//^^^'^ and one gets, 
at 95% confidence level, 

log (^^^l^j < 0.52. (127) 

The marginalised probability distribution for the reheating parameter is peaked 
around In R ~ —22. Its behaviour at large values directly comes from the prior (11261) 
which is a function of pstrg- As will be discussed in the following, pstrg is related 
to the energy scale of inflation which is bounded from above by the amplitude of 
the cosmological perturbations. Here again, the differences between the marginalised 
probability and the mean likelihood trace the correlations induced by the prior 
hypothesis between these two parameters, namely that reheating occurs after inflation. 
The lower tail of the distribution is driven by the data and falls off till the prior lower 
bound is saturated at In R = —46 (nucleosynthesis limit). This posterior is very similar 
to the one derived for the small field models in reference [M] where it has been shown 
that the CMB data were disfavouring a low energy scale reheating for these models. In 
the present case, the In R posterior of figure [8] does not fall-off to zero in its lower part 
but still yields a limit 

lnR>-38, (128) 

at a 95% confidence level. 

The physical interpretation of this limit comes from the influence of R on N^, 
the e-fold at which an observable perturbation mode with wavenumber k^, crossed the 
Hubble radius during inflation. Indeed, one has a^,/ao = a^:/astrg x astrg/c^O) with, by 
definition, a^^/ao — K/{aoK,^/^V^^'^), where we have expressed the Hubble parameter 
in terms of the potential. The quantity a=„/astrg can be expressed in terms of A^,,, and 
Ostrg/flo = (1 + -Zstrg)"^ is glvcu by equations (11191) . (11201) and (I12ip . expressing pstrg in 
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terms of the potential evaluated at 0strg- One gets 



iV, ~ 58 - In 
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As shown in section [71 since the WMAP data prefer a slightly red-tilted spectral index 
in the absence of tensor modes, it is not surprising that, for a given model of inflation, 
some values of iV*, and therefore R, end up being favoured. This interpretation can 
be further explored by using the slow-roll approximation. As already mentioned, the 
KKLMMT model generically leads to observable ei values much smaller than €2- Let us 
first determine determine iV^, such that €2(0*) = e2obs; where (j)^ = (p{N^) and e2obs — 0-05 
is the preferred observed value (see figure [7]). Using the classical trajectory given by 
equation (J5T|) . one finds that 
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and using the expression of the second horizon flow parameter (!50l) . one arrives at 
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Combining these two last equations, one obtains 
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As already discussed, if the slow-roll is violated before brane annihilation, we can use 
the slow-roll trajectory with the approximation 0end — 4'e2 — [40/(/«/U^)]^''^, from which 
one deduces that 

^ . (133) 
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Inserting this value into formula (11291) with k^/ao ~ 0.05 Mpc^^ and neglecting the 
potential term gives In R ~ —28. This is the preferred value of In R that would give an 
observable spectral index at /c* compatible with the favoured value and in agreement 
with the shape of the marginalised posterior of figure [Hi Let us notice that the two-sigma 
limit (11281) is the best that can be extracted from the data. Indeed, the three-sigma lower 
limit matches the nucleosynthesis prior bound. This differs from the results obtained 
in reference [51] for the small field models and is due to the marginalised posterior 
which does not vanish for low value of In R. Moreover, the mean likelihood appears to 
increase again in that region. This effect comes from correlations between parameters 
and suggests that some, but quite fine-tuned, models provide a good fit to the data even 
for In R small. To understand this property, we have represented in figure [HI the two- 
dimensional marginalised posterior and mean likelihood in the plane [In R, log(0strg/ A*)] 
(left panel). One recovers the favoured value for In R ~ —22 as the central vertical 
broad region (one-sigma contour), and as shown above, this domain is associated with 
the models in which slow-roll inflation ends before the branes collide. The two-sigma 
contour extends towards the left side of the plot but only for 1 < log (0strg//^) < 1-5. 
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Figure 9. The left panel shows the one- and two-sigma confidence intervals of 
the two-dimensional marginalised posterior as a function of the parameters In R and 
log{(f)strg/ fJ-)- The shading traces the corresponding two-dimensional mean likelihood. 
The right panel is a contour plot of In R as a function of log{^/Kfi) and log((/)strg//^) 
obtained from the slow-roll approximation by assuming that 62 — e2obs — 0.05. The 
dashed line marks the limit (/)strg = (t'e2- This plot shows that, although (/)strg > M 
is generically disfavoured (see figure [5]), it is still possible to obtain a good fit to the 
data provided some fine tuning is performed between (pstrg/fJ^, /^/'tt-pi and In R: a low 
energy reheating allows for some models where brane annihilation occurs before slow- 
roll violation to have red-tilted power spectra. 



This tail is associated with the models in which brane annihilation occurs before the 
slow-roll conditions are violated. To see why such models are not disfavoured we can 
still use equations (11291) and (11321) but, now, with 0end = 0strg to get 



The In R contour plot coming from this function is represented on the right panel of 
figure [9] for different values of 0strg and ^/rupi. One sees that it is indeed possible to 
get e2obs — 0.05 for 0strg > 4>e2 provided In R ~ —40, in agreement with the plot. The 
physical interpretation is that a low energy reheating allows the observable window 
of wavenumbers to be shifted towards the end of brane evolution in such a way that, 
although brane annihilation occurs before violation of the slow-roll conditions, it is 
possible to get not so small values for 62- Of course, this requires some amount of fine 
tuning between In R, n/mpi and 0strg; which is statistically penalised and explains why 
these models are out of the one-sigma contour. 




(134) 
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8. 6. Derived primordial parameters 

As described in section Ell the Markov chains have been performed on the power spectra 
amphtude and its posterior probabihty distribution is plotted in figure El Since V* 
and M/rripi are in direct relation (see section [8l2|) . it is possible to derive the marginalised 
probability distribution of M/ rripi from the one of . The same considerations apply to 
the volume ratio v which is given by equation flMl) : using importance sampling, one can 
extract its posterior distribution from the one of /i/mpi and M/m^i [7T]. Both probability 
distributions are plotted in figure [10] and seem to favour some peculiar values of M/ mpi 
and V. Once again, the discrepancies between the posteriors and the mean likelihoods 
come from the correlations with other parameters and may hide some effects coming 
from the priors. The right panel of figure [10] represents the two-dimensional probability 
distribution (point density), as well as its one- and two-sigma confidence level regions, 
obtained without marginalising over the parameter \og{^/K,^). As can be seen on these 
plots, the narrow highly probable regions trace the strong correlations between M/ rripi 
and fi/nipi, and between v and fi/nipi. Remembering that the lower limit on the /x/mpi 
parameter comes from our numerically convenient prior y/nfi > 10~^, one immediately 
sees that this choice has a direct infiuence on the upper limit, and respectively the lower 
limit, of the log(47r^f) and log{^/K,M) probability distributions. If we had chosen a 
smaller limit for fi/rripi, the v posterior would have been shifted towards larger values, 
whereas the M/nipi one towards the smaller. As a result, we conclude that there is no 
upper, respectively lower, constraints on v and M/nipi. 

The same correlations are at work for the other tails of the probability distribution 
associated with v and M/rripi. However, this time, we are probing the highest allowed 
values for fi/nipi which are disfavoured due to the amount of fine-tuning required for a 
successful infiation nearby the throat edge (see previous section). As a result, we obtain 
the 95% confidence levels 



Let us stress that these bounds come from both the data and the requirement that 
infiation proceeds inside the throat, as imposed by the self-consistency of the model 
used. 

To understand how the data lead to the correlations observed between the 
parameters fi/rupi, M/mp^ and v, we have plotted in figure [TT] the one- and two-sigma 
contours (solid curves) of the two-dimensional posteriors associated with various pairs of 
primordial parameters. The corresponding two-dimensional mean likelihood is traced by 
the shaded areas. The highly probable region spanning the plane [\og{y/KM) , \og{y/Kfi)] 
directly comes from the constraints on the power spectrum amplitude and spectral index. 
This can be shown by using the slow-roll approximations. Indeed, if one assumes that 
the slow-roll conditions are violated before brane annihilation, namely (pstrg < <Pe2i then 
the WMAP normalisation through equations f[68|) and ([69]) leads to 



log (v^M) < -2.9 



log {A'k'^v) > -8.5 . 



(135) 



log (v^M) = log 
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Figure 10. Marginalised posterior probability distributions for the M/mp, and 
V parameters (solid curves) and their associated mean likelihoods (dotted curves). 
The right panel shows the corresponding one- and two-sigma contour of the two- 
dimensional posteriors obtained without marginalising over \og{^/Kfj,). The two- 
dimensional probability is proportional to the point density while the colormap traces 
correlations with a third parameter. 



+ ^log(v^/i). (136) 

This is a straight hne with slope equal to 1/3, exactly as observed for the high probable 
regions in figure [TTl Using a fiducial value for the quadrupole and A^,,, = 50, the offset 
is around —2.6. Consequently, for log (v^/^) — ~3 one gets log (a/kM) ~ —3.6, as 
confirmed by the plot. The effect of the spectral index is indirect and, as previously 
discussed, tends to favour the models in which (p^trg < 4>e2- This suggests that the less 
probable and wider shaded region in the plane \[og{^/KM),\og{^/R,fJ,)] corresponds to 
the cases where 0strg > • Again, the slow-roll approximation in this limit gives the 
relation 



Since v is also a function of M/mpi and fi/rripi, this is a three parameters relation, 
or a surface in the plane [log{y/K,M),\og{y/Kn)]. We can, however, derive the slope 
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Figure 11. One- and two-sigma contours (solid curves) of tiie two-dimensional 
marginalised probability distribution for various pairs of the primordial parameters. 
The two-dimensional mean likelihood is traced by the shading intensity. The bottom 
panels show the two-dimensional posterior (point density) and its correlation with a 
third parameter (colorbar). 



associated with the lower boundary of the blurred region (beyond which there is no 
acceptable models). This boundary is reached for the largest 0strg compatible with the 
throat size, i.e. for 

in,strg — max(0uv) — f^pi/ 'n/Stt [sec equation f l35l) ] . This limit is 
reached when the ratio of throat to bulk volume is maximal. Using the expression of 
</'in,strg given by equation (IMl) and solving for 0strg leads to 

0strg 
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This expression can now be inserted into into equation fll37p and one gets 




5 , 64 24A^T 



(139) 



12 {^f,f ^ 



For not too large values of a/k/x (in practice a/k/^ ^ 0.1), the term proportional to A^^ 
can be neglected and one obtains 



This is a straight line in the plane [log{y/K,M) ,\og{y/K,fi)] the slope of which is +2: this 
matches with the observed lower boundary of the blurred region in figure [Til The offset 
is around —2.1 giving \og{^/KM) ~ —4.1 for log (a/k/i) ~ —1, again in agreement with 
the figure. One might also notice that the boundary curve in figure fTTl slightly bends over 
for small values of fi/rripi as the effects of the terms proportional to N^, start to appear. 
The previous interpretation is confirmed by the left bottom panel of figure [TTl The 
point density clearly decreases in the domain of low likelihood which appears blurred in 
the upper left panel. As traced by the colormap, this region indeed corresponds to large 
value of 0strg/Ai- Let us notice that, although this region is out of the 95% confidence 
contour, it remains inside the three-sigma one. 

The previous discussion also applies to the right plots of figure [TOl The above- 
described relations between M/nipi, /x/mpi and 0strg are directly converted into 
correlations between v and the other parameters through equation fl2^ . As a result, the 
main degeneracy seen in the plane [log(A/K/i), log(47r^f )] is also a consequence of both 
the power spectra normalisation and spectral index constraints. The bottom right panel 
illustrates the degeneracy between M/rripi and yu/mpi coming from the data and differs 
from the one associated with equation (124]) : this difference explains the two-sigma lower 
limit on log(47r^f). Let us elaborate on this point. The limit on v can be understood 
from the slow- roll result even if, at first sight, there is a mismatch between the constraint 
obtained above and equation flllll) . In fact, equation (11051) implies 



This relation is consistent with the top right panel in figure ITOl Let us notice that, in this 
context, it is relevant to use equation (11051) . which is derived under the assumption that 
inflation stops by violation of the slow-roll conditions, since one observes in figure [lO] 
that the two-sigma contour corresponds to "small" values of 0strg- Moreover, if one 
inserts the limit fl89l) . namely y/^fi < 2, in the above relation then the constraint (11351) 
is reproduced. This constraint is different from equation (lllip because it has a 
different origin. The limit (II 111) assumes that the maximum allowed contribution of 
tensor modes to the observed CMB data can indeed be generated during KKLMMT 



log (a/kM) ~ log 45 




+ 21og(v^/i) . 



(140) 



log (^An'^v) — — - log (v^/i) — 9 . 



(141) 
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inflation. However, the requirement that the brane motion proceeds inside the throat 
strongly limits the generation of tensor modes and leads to the stronger bound of 
equation f llSSp . Therefore, the two limits are consistent and the stronger bound is 
given by equation fllSSp . 

Finally, we have plotted in the right hand panels of figure [TT] the two-dimensional 
probability distribution and mean likelihood in the plane [log{y/K fi), logi^cpstrg/ fJ^)], as 
well as the effect of (pyy. These plots exhibit the volume effects associated with the 
large values of the parameter /i/mpi. Clearly, the allowed range of both 0strg/At and 0uv 
are all the more reduced as /i/mpi increases, thereby decreasing the statistical weight 
of these domains in the marginalised posterior of yu/ mpi (see figure El) . The two-sigma 
contours are found to follow the high likelihood region, and as before, it corresponds 
to models in which brane annihilation occurs after violation of the slow-roll conditions 
(0strg < (pt2)- Indeed, the edge between the dark and the blurred regions is just given 
by the equation 0strg = 4>t2- Using equation (16T!) . which gives when a/k/x is small, 
one gets 



This is a straight line in the plane [log(A/fi;yu), log(0strg//^)] with a —1/3 slope, as 
observed in the figure. The offset is close to 0.3 implying that for \og{^/Kfi) = —3 
one would have log(0strg/Ai) ~ 1-3. There is a slight difference since the contour seems 
to intercept the vertical axis at log (0strg//^) ~ 1-5. Such a difference may be due 
to the fact that (/)strg must significantly deviate from fo^' the effect to be visible. 
Moreover, the reheating effects already mentioned are not considered in the derivation of 
equation f ll42p . Equation (11420 renders explicit the //-prior dependence previously noted 
on the upper two-sigma limit associated with the marginalised probability distribution 
of 0strg//^ (see figure (HI). The —1/3 slope explains why the prior independent limit of 
equation (I127P is more conveniently express in terms of (j)/ fj?^^. Finally, the upper edge 
of the blurred region is, as before, given by the condition 0in,strg = niax(0uv) (see lower 
right panel of figure [H]) . Using equation (I138p , one finds 



Again, one can check that the slope and the offset are consistent with what is observed 
in figure [TTl Beyond that limit, infiation would not take place in the throat and those 
models are rejected. 

8. 7. Fundamental parameters 

Our choice of MCMC parameters is minimal as far as observable effects are concerned. 
The cosmological consequences of the KKLMMT model are described by the set of four 
parameters (M, /i, 0strg, 0uv) whereas, as discussed in section 12.51 the underlying string 
theory model involves five parameters: {gs,a',M,v,M'). As a result, it is impossible 
without additional assumption to extract more information for the string parameters. 
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Figure 12. Marginalised probability distributions on the reheating energy scale 
derived from the WMAP data under the assumption that the reheating proceeds with a 
constant equation of state: P = WrehP- For an extreme equation of state Wreh ^ —1/3, 
one finds a weak (non-trivial) lower bound > 20 GeV at two-sigma level. For all 
the other cases, the reheating can occur at any energy scale higher than nucleosynthesis 
(from the CMB point of view). 



As found in the previous section, the CMB data, however, aUow some constraints to be 
derived for the M/rUpi and v parameters. Their probabihty distributions do not assume 
anything on the value of Qs, a' or H (apart from consistency requirements) and are also 
robust against the reheating as long as it can be described by our phenomenological 
model. If one wants to go further on the fundamental theory parameters, additional 
assumptions have to be made. In the following, such a step is performed first on 
the reheating model by assuming that it proceeds with a constant equation of state 
P = Wj-ciip. Then, a similar approach is adopted to derive some posterior probability 
distributions on the remaining string parameters by assuming that the vahic of a' is 
known. As a working example, a' = lO'^mpj^ and a' = 10^m~^ are considered. 

For a constant equation of state parameter during reheating, p oc a-3(i+«;reh) g^^fj 
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equation fll2ip can be further simplified into 

From the MCMC analysis performed in the previous section, samples on R and pstrg can 
be used to extract by importance sampling the probability distribution associated with 
Preh (the prior choice remains that the reheating occurs before nucleosynthesis). The 
resulting posteriors are plotted in figure [T5] for four values of Wreh spanning the range 
allowed by the dominant and strong energy conditions in General Relativity. For an 
extreme equation of state with w^eh ^ —1/3, i-e. which is on the verge of an accelerated 
expansion, one has the 95% confidence limit 

p;4'>20GeV. (145) 

Although this kind of models are certainly already ruled out by particle physics 
experiments, let us notice that this limit comes from CMB data only. As can be seen 
in figure [121 ciU the other cases associated with more reasonable values of the equation 
of state parameter are not constrained. 

A similar method has been applied to the parameters Qs and by assuming 
that the value of a' is known. In figure [13], we have plotted the two-dimensional 
marginalised probability distributions, as well as their one- and two-sigma contours, 
in the plane (logA/", logo's) and for the two fiducial values of a'. As can be seen on these 
plots (right panels), there is basically no constraint coming from the data for the case 
a' = 10^m~^. Indeed, on one hand, the two-dimensional posterior is limited by the 
consistency conditions Af > 1 and (7s < 1 to satisfy flux quantisation and perturbative 
treatment. On the other hand, the colorbar shows that the two-sigma contour simply 
traces the high values of v which have been shown in section [8l6] to be essentially related 
to our lower prior on /i/mp,. The case a' = 10^m~^ is slightly different. As can be 
seen in the lower right panel of figure [131 the two-sigma contour is detached from the 
axes and therefore the data give slightly more information than the prior consistency 
limit. Again, the right end of the contour is not physical since it comes only from our 
fi/nipi prior, as before. However, the observed degeneracy between log TV and logo's is a 
non-trivial information. It means that for a given flux number, a value of Qs is favoured 
by the data. We have not represented any one- dimensional posterior for Qs or Af since 
there is typically no more information than the above-mentioned degeneracy. To do so, 
another assumption would have to be made, such as some prior knowledge on the flux 
number or on the coupling constant. On the left panels of figure [131 ^6 have plotted 
the corresponding two-dimensional posteriors in the plane {ln(47r5fsA/'/f ), ln[f /(47r(7s)^]} 
to allow a comparison with the slow-roll results of section [6] (see also figure [5]). As 
expected, all the models lie in the range predicted by the slow-roll analysis, up to the 
prior limit fixed by the numerically convenient lower limit on fi/rripi. The associated Qs 
values are traced by the colormap. 
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Figure 13. Two-dimensional posteriors of the string parameters and their one- and 
two-sigma contours obtained by importance sampUng for two fiducial values of a': 
a' = lO'^mpi^ and a' = lO^mpj'^. Correlations with a third parameter are traced by 
the colormap. The CMB data give some non-trivial constraints only for large values 
of a' . This appears as the degeneracy observed between J\f and gs in the bottom right 
panel. The left panels show that the models lie in the regions predicted by the slow-roll 
analysis in section [6l 

9. Conclusions 

To conclude, let us summarise and discuss the main results obtained in this article. Our 
goal was to compare a typical model of string-inspired brane inflation to the CMB data. 
In this scenario, the inflaton fleld is interpreted as the distance between two branes 
moving in six dimensions along a warped throat. The end of inflation either occurs 
by violation of the slow-roll conditions or by tachyonic instability at brane annihilation. 
Our MCMC analysis was carried out by using an exact numerical approach in which the 
only approximation used is the linear theory of cosmological perturbations. Moreover, 
various effects have been considered as, for instance, the presence of a DBI kinetic term 
or the possible quantum effects. The MCMC exploration has also been performed for 
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arbitrary values of gs and a'. 

From this analysis, we have obtained the following results. Firstly, the WMAP 
data favour a scenario where inflation ends by violation of the slow-roll conditions 
before brane annihilation rather than by the tachyonic instability brought about by 
the annihilation. In other words, 0strg cannot be too close to the edge of the throat and 
one finds 

-1/ 



log < 0.52 , (146) 



^2/3 

at 95% confidence level. This constraint originates from the fact that ending inflation 
by instability pushes the spectral index towards one while preserving a low level of 
gravitational waves, a situation which is disfavoured by the data. 

Secondly, one has obtained a limit on v, the volume ratio of the five-dimensional 
sub-manifold forming the basis of the six- dimensional conifold to the volume of the 
five-sphere. This limit reads 

logt;>-10, (147) 

at 95% confidence level. This constraint comes both from the data and a requirement 
of self-consistency of the model, namely that inflation proceeds inside the throat. 

Thirdly, we find that the reheating period is slightly constrained. Although we have 
not considered the detailed process of brane annihilation, the total number of e-folds 
for which the universe reheats may change the part of the inflaton potential probed by 
the observations. Combined with the power spectra generated during brane inflation, 
the WMAP data provide the 95% confidence limit 

In R > -38 . (148) 

In the case where the reheating proceeds with a constant equation of state parameter, 
we find that if t^reh ^ — 1/3 then pl^^ > 20GeV (at two-sigma). This limit is certainly 
already ruled out by particle physics experiments, but it is worth recalling that it has 
been obtained from the CMB data only. With more accurate data, it will certainly be 
possible to improve this bound, and hopefully for other values of the equation of state 
parameter. 

Fourthly, on the theoretical side, we have obtained approximate solutions in the 
case where the DBI kinetic term does not reduce to the standard one. The regime 
where the initial conditions are such that the quantum effects are important has also 
been considered in detail, see Appendix A[ In this case, using a perturbative treatment. 



we have computed how the trajectory (with or without volume effects) can deviate from 
the classical motion. We have also shown that this approximative scheme breaks down 
when the brane approaches the bottom of the throat. In this regime, only a numerical 
integration of the Langevin equation could allow us to go further and, for instance, to 
see whether the field really starts to climb the potential as indicated by the calculation 
presented in the appendix. Clearly, this discussion is relevant in order to know whether 
a regime of eternal inflation can be established in this model. 
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Finally, let us discuss how the present work could be improved. Recently, various 
works have been devoted to the type of model studied here. In particular, special 
attention has been paid to the general set up and exactly calculable corrections to the 
potential (fT3!l have been proposed. It was also show that, in some situations, these 
corrections play a crucial role [2ll ES, EEl ETJ EH]- The next step would therefore 
be to include them in our analysis. Since they introduce new parameters into the 
problem, and given the weakness of the constraints obtained here, this would probably 
be meaningful only at the expense of including other more accurate data sets to break 
the degeneracies. Time dependent KS-like compactifications have also been discussed 
in the literature to provide a late-time acceleration of the universe [9^ 195] . These 
models may provide an alternative to the standard ACDM universe we have considered 
at low energy. Concerning the approach of the present paper, more precise data would 
certainly help to disentangle the correlations between J\f and Qs, and might therefore 
allow to use astrophysical data to constrain String Theory, at least under the strong 
theoretical prejudice that the relevant model of inflation is indeed the one described 
here. 
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Appendix A. Stochastic inflation 

Appendix A. 1. Regime of quantum fluctuations within the throat 

We previously stated (section [231) that all e- folds of infiation should occur within a single 
throat and hence < always, being the field value at the edge of the throat. In 
section [3l we further calculated the field value 0fiuct above which quantum fiuctuations 
have an impact on the classical trajectory. It seems reasonable to ask under which 
conditions, if any, the field value 0fluct is located within the throat and hence stochastic 
infiation effects can, at least in principle, affect the field evolution. Let us try to derive 
the parameter restriction resulting from 

0fluct < 0UV • 

j| http://www.idris.fr 
% http://www.planck.fr 
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Figure Al. The panels illustrate the parameter regions allowing the stochastic regime 
to take place inside the throat, for a'm'^^ = 1000 (left panel) and a'mp, = 10 (right 
panel). All curves are the same as in figure [5l but in addition a solid green curve 
illustrates the condition (/ifluct = (t'vv- Stochastic inflation may occur inside the throat 
only for the parameter values lying above the solid green curve curve, and below the 
volume constraint (dashed green curve with ticks down). 



Again, it will be important to decide which normalisation is used to express 0fluct given 
by equation (HGl) . With from equation (flQl) . consider first the normalisation (l68ll 
valid for 0end = <Pe2- this case, the rescaled parameters (x, v) of equation (1731) are 
again a useful set to express flA.l|) in a simple form, 

Int; > 41n J"- Inx, (A. 2) 

where the constant can be expressed as 

^ 3I/102-1/5 /45^^£^j (6Ar, + 10)-^/' (Wm^y/^ . (A.3) 

This is a condition of the same type as the one given by equation (19^ (from 0in,e2 < 0uv)' 

but more restrictive. The parameter region in the [Inx, \nv) plane for which stochastic 

effects can occur within the throat is therefore smaller than the one obtained from 

equation (p2l) (see figure ETj) . 

Now let us investigate the case of equations (!99l) and ( llOOp where the WMAP 

normalisation is done using 0end = 0strg ^ • Inserting this normalisation into 

equation (H6l) for ^Auct? we find the following inequality in terms of (Inx, In -y): 

5 12 
In^ < — lnx + - lnw + (A.4) 
12 3 3 ^ ' 



where 



g ^ 3I/IO2-I/5 (^c.'^2 J 1/2 9!l^^ . (A.5) 

This is another curve in the plane (Inx, Inv), analogous to fllOip though slightly simpler. 
Again, equation (1A.4P constrains the possible parameter space outside the contour 
0^2 = 0strg in the same manner, but somewhat tighter than equation (llUlj) . Note 
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that inside and outside of the contour (j)^^ = 0strg, the requirements 0in,£2 < 0uv 
0fluct < 4>vv consistent with each other. Both are represented in figure IA1[ We have 
therefore shown that the condition ( lA.ll) can be met for sufficiently generic parameter 
combinations, and hence stochastic effects may occur inside the throat. 



Appendix A. 2. Detailed analysis of the stochastic regime 

In the following, we apply the treatment established in reference j60] to the case of 
the brane inflationary potential f|T3l) . The underlying approach of stochastic inflation 
consists in describing the evolution of a coarse-grained fleld ip corresponding to the 
original scalar fleld averaged over e.g. Hubble patch. The quantum fluctuations 
correspond in this picture to stochastic noise due to small-scale Fourier modes. 
Therefore, the dynamics of the coarse-grained fleld are controlled by a Langevin 
equation, whose solution enables one to calculate the probability density function of 
the fleld (p. In the slow-roll approximation, this coarse-grained fleld obeys the equation 

where the Hubble parameter H is entirely described by the coarse-grained fleld as well, 
i.e. one has (in the slow-roll approximation) 

^ ^Vi^if). (A.7) 

In equation (1A.6I) . ^(t) is the noise fleld and its mean and two-point correlation function 
satisfy 

(^w) = o, mm) = Kt-t'). (A.8) 

where 5 is the Dirac delta function. The approach in jBOj has been to solve this system 
of equations using a perturbative expansion in the noise, that is, one writes ip as 

Vit) = Mt) + Sviit) + SMt) + • • • , (A.9) 

where (/)ci(t) is the solution of the Langevin equation without the noise, and S{pi{t) is 
linear in ^{t), Sip2{t) quadratic and so on. To avoid notational confusion, the classical 
fleld value, which is denoted by in the rest of this article, is called 0ci in this Appendix. 
As discussed in reference [60], there are, up to second order, two resulting probability 
distributions for ip. The flrst of these, Pc{p, t), is the probability of the stochastic process 
to assume a given fleld value at a given time in a single coarse-grained domain. The 
second distribution, P^{p,t), takes into account volume effects, namely the dependence 
of the volume size on the fleld value within this domain. These two distributions are 
given by 



P,{p,t)= {6{p-ipm= exp 

V27r {5ipi) 



2 {6^1) 



(A.IO) 
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and 



(%-y.[e])exp[3/drg(y.[e])]) 
(exp [3fdrH{^[^])]) 



: exp 



(A.ll) 



2 {SvD 

(pel + {^^2) and where / and J are 



where the mean field value is given by ((/?) 
two continuous vectors [60]. The infiuence of stochastic effects therefore requires three 
quantities to be determined for any given inflationary potential, namely {Sipf), {Sip2), 
and I^J. These quantities have been calculated in reference [60] as functions of the 
Hubble parameter (in the slow-roll approximation) and read 

2 ^ih:„ / TT\ 3 



K 

2 

H" 



H' 
2^ 



<t>cl 



H 



2H 



7 m + 



H' 



Hi 



I 2 



H' 



12 



(A. 12) 



(A.13) 



where, in the present context, a prime denotes a derivative with respect to (p. In the 
slow-roll approximation, the Hubble parameter can be directly replaced by a/ kV/S and 
the correction due to volume effects can be expressed as 



3J^J 



12H' 



ml 



127r 



H' mi. 



(A. 14) 



In the stochastic inflation regime, we are dealing with very large fleld values ^ /i 
and one can use the expansion f[T^ of the potential. Equation flA.12p then yields 

-1 

12" 



(A.15) 



Since the brane motion in the throat occurs only for 0ci > (see section Hj), one gets in 
the regime 0ci ^ 0in 

^'''^ . (A.16) 




7687r2 



Since 0ci//^ is decreasing, the variance is, as expected, increasing as inflation proceeds. 

Let us now determine the correction to the mean value. As can be seen from 
( lA.lSp . no additional integration is necessary to calculate {S(f2), and after some algebra 
one flnds 



0cl 



0cl 



rii _ 1 
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X 



1927r2 



1 
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(^ci) _ 


-1/2 ^ 










1 5/2 



t^ 
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5/2 ^ 
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(^ci) _ 
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(A.17) 



Finally, in order to evaluate the volume effects characterised by the term 31^ J, one 
needs to calculate the integral in equation (1A.14I) . For the potential (fT4l) . this can still 
be solved exactly. One obtains 



1287r2 



^cl 



1 - 



5 



S 



1 



0cl 



^1 



(A.18) 



where the function S is defined by 

1 



S{x) 



+ 



383 
105 In 



1 (-279 + 326a;^ - 200x^ + 48x^^) 



a; 



(A.19) 



The corresponding trajectory is plotted in figure |A2]( left panel). The solid line represents 
the classical trajectory without the quantum effects. Since we are using 0ci as a clock, 
this line is just the straight line = x". The dotted line gives (y?) = 0ci + {^^2), i-e. the 
mean value of the field in the situation where the volume effects are ignored. The two 
dashed lines above and below this line are 0ci + {^^2) ± ^/ {^^D- As discussed at length 
in reference [60] , the quantum effects accumulate and manifest themselves only once the 
field is approaching the end of inflation. As clearly shown in figure IA21 the quantum 
effects are such that the field rolls down the potential more quickly than its classical 
counterpart. However, the story is different when volume effects are taken into account. 
The corresponding trajectory is represented by the dotted-dashed hne, the two dashed 
lines above and below being the limits of the one sigma interval. In this case, the field is 
first slowed down and starts moving back to even climb up the potential. Then, at some 
point, it stops and rolls down again the potential. However, as we will see in the next 
subsection, this last part of the evolution can not be tracked within our approximation 
since it is outside its domain of validity (represented by the shaded regions). 

In figure IA2I (right panel), we have represented the corresponding distributions 
without and with volume effects as given by equations (lA.lOl) and flA.lip respectively. 
The vertical dotted lines corresponds to the classical values of the field at three successive 
times (denoted "1", "2" and "3"), namely 0^fV/^ = 300, = 200 and (p'^'^/fx = 150. 

The solid lines are the three distributions Pc corresponding to the three classical field 
values while the dashed lines are the three Pv's. This plot essentially confirms the above 
analysis. In particular, one notices that the maximum of Pc is ahead the classical vacuum 
expectation value as discussed before. On the contrary, the volume effects "retain" the 
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400 



Figure A2. Left panel: Evolution of the inflaton field when the quantum effects (with 
or without volume effects) are taken into account. The underlying model is such that 
"classical" inflation stops by violation of the slow-roll conditions. Hence, the WMAP 
normalisation is calculated with the help of equation ([68]) . With w = 1 and N^, — 50, 
this leads to (pauct/f^ ~ 793.5. The initial condition is chosen to be (pin/ fJ- — 700. The 
solid line represents the classical trajectory, the dotted line the mean value without 
the volume effects and the dotted-dashed line the mean value with the volume effects. 
Right panel: Evolution of the probability distribution Pd^fi) (solid lines) and Pv{<p) 
(dashed lines). As in figure \X2\ the case displayed here corresponds to u = 1 and 
iV* = 50, leading to (f>finct/l^ ^ 793.5. The initial shape of the probability density 
function has been chosen to be 6{(j) — (pin) with (pin/ ^J- — 700. The vertical dashed lines 
represent the location of the field at successive points of the classical trajectory (the 
field moves "from right to left" ) . 



field and this one can even change the direction of its motion (from right to left to left to 
right). Finally, the two distributions strongly spread out while inflation is proceeding. 



Appendix A. 3. Reliability of the perturhative treatment 

In this section, we discuss in which range of field values {ip) G [0ci ~ |Av^min(0ci)|, 0ci + 
A<^max(0ci)] the previous analysis is valid. Indeed, as discussed above, in order to 
estimate the mean value of the field with or without the volume effects, a perturbative 
expansion in the noise has been used which may break down at some point. It has 
been shown in [61] that Ay^min and A(y9max can be found from requesting that the two 
conditions 



max 

n 



max 

max 
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TTlll 



X 



H. 



3/2 



Alp 



(A.20) 
(A.21) 



are simultaneously satisfied. The strongest bounds on A(/?min and Ayj^ax that follow 
from (IA.20p and (IA.21I) then give the reliability interval for the distributions Pc{'p>it) 
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and P^{!f,t). In the present straightforward calculation leads to 

A^mm^_2 0cl A^max _ 2 0cl (A 22) 

The resulting reliability region is plotted in figure IA2I and corresponds to the purple 
hatched region. We see that the approximation breaks down more quickly for Py[(p,t) 
than for Pc{(p,t). In particular, the conclusion that, due to volume effects, the field 
reverses its motion comes from a regime which is at the border of the reliability interval. 
We conclude that, in order to have access to this phenomenon and, in particular, 
to discuss eternal inflation in this model, a more powerful method is needed, as, for 
instance, a full numerical integration of the Langevin equation (1A.6p . 



Appendix B. Brane inflation from String Theory 

Appendix B.l. The background action of type IIB Superstring Theory 

We begin by considering the background action only, and to that end place ourselves 
in the context of ten-dimensional type IIB Superstring Theory. In the string frame, the 
effective action reads (see e.g. references [i9 | \57\ [96j) 



, _ 1 

^tot — 



2k 



/? + 4(V$f -l(di?(2))^ 



+ C(4) A dfi(2) A dqs) , (B.l) 
where k stands for 

k=^-{2nya'\ (B.2) 
The fleld strength denoted by F(5) in equation flB.lD is given by 

i^(5) = dC(4) + I [Bi2) A dC(2) - C(2) A dS(2)] . (B.3) 

The brackets (z) indicate the number of indices of the corresponding fleld. The flelds 
in this theory are therefore the scalars $ (the dilaton) and C(o) (the axion), (ten- 
dimensional gravity) as well as the antisymmetric flelds 5(2) (Kalb-Ramond), C(2), and 
C(4). Note that the totally antisymmetric e tensor refers to curved space, where we 
deflne 

6(10) = ±1 for pair /odd permutation of the indices, (B.4) 
^(10) =0 if two or more indices are equal, (B.5) 



e('°) = -6(io). (B.6) 
9 

The flrst line of ( IB. 11) describes the Neveu-Schwarz-Neveu-Schwarz (NS-NS) sector, 
while the second line results from the Ramond-Ramond (RR) sector. The third line 
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is the Chern- Simons (C-S) couphng between the sources of C(4) (that is, RR charged 
objects such as D3 branes) and background fluxes. From the definition flB.3|) of -F(5) it 
is clear that if (d-B(2)) and (dC(2)) have non-zero background values, they can produce 
a non-zero F(5) field strength even in the absence of locahsed C(4) sources. Turning 
this argument around, one can also interpret such fluxes as intrinsic "background" C(4) 
sources. 

Appendix B.2. The equations of motion 

From equation ( IB.ll) . it is straightforward to derive the equations of motion for these 
fields. The resulting system of equations has been studied in the references [ISlBQllSm iST] 
and we now present the solutions found there. Note that in addition to finding 
simultaneous solutions for all equations of motion, one has to impose the condition 
of self-duahty of the 5- form field strength, F(5) = ^F(5), onto these solutions. 

Appendix B.2.1. Ansatz for the metric. Let us focus on the Einstein equations 
first, which follow from ( IB.ip by variation with respect to the ten-dimensional 
metric (^mnEI- In this section, the ten-dimensional set of coordinates we use is 
x\ r, ip, 6*1, 01, 02, 02) , ^ = 1, 2, 3. A super-gravity motivated ansatz for the metric 
reads 

ds' = v.udx^dx" + ^/h(^dsl , (B.7) 

^/h{r) 

which involves a natural splitting into a four-dimensional extended space-time and the 
six compactified extra dimensions. Note that in (IB. 70 . the four-dimensional metric has 
been set to Minkowski r/^,^ = diag(— 1, +1, +1, +1), multiplied by a factor depending on 
the fifth dimension r. In principle one can also choose e.g. a FLRW universe for the 
four- dimensional section of the metric by replacing rj^^, with g^^, = diag (—1, a^, a^, a^). 
For the six- dimensional section, the choice of interest here is 

ds^ = dr^ + r^dslii, (B.8) 

where ds|.^ ^ is the metric on the Einstein space Ti i [17]. Note that the choice of 
geometry and compactification for the extra dimensions is an input parameter for the 
resulting inflationary model. In particular, one can enforce "warping" on the conifold 
Ti l by giving non-vanishing background values to (di?(2)) and (dC(2)). The main virtue 
of these fluxes is that in general they stabilise the scalar fields corresponding to the 
complex structure moduli of the extra dimensions. The warped conifold geometry was 
studied in reference [50] and has become known as the Klebanov-Strassler (KS) throat. 
Interest in the KS geometry is due to the fact that it represents an exactly solvable 
example of field strengths background fluxes. Nevertheless, this scenario can be viewed 

+ Small Latin indices m,n,. . . run from to 9 (with the exception of i = 1,2,3), while small Greek 
indices /i, t^, . . . run from to 3. Capital Latin indices M, N, . . . run over the extra dimensions 
coordinates only, i.e. from 4 to 9. 
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as generic because "warped throats" appear in many flux compactifactions invoked to 
stabilise the moduli. The metric of the Ti i space reads 

= dr^ + — {d-ip + cos 9i d0i + cos 62 d</)2) 
9 

^2 



1=1 

whereas the one on its deformed counterpart is 

_ _ 1 ,0 r2 

'6, warped " ^(-^^ ' g 

^2 



dSg, warped = + —K{r) {dip + COS ^1 d0i + COS ^2 ' 

„2 



o 

+ i (r^ + 6z^) {del + sin^ ^2 d^j) , (B.IO) 
where the warping function K{r) is given by 

7-2 _|_ 0^2 

A-M-;^. (B.ll) 

The quantity z 7^ (in general complex) is the warping number introduced in the 
definition of the deformed conifold (see below). 

We begin our investigation with the simple conifold (IB.Op . Using our coordinates 
(r, ip, 9i, (pi, 6*2, 4>2) for the six- dimensional section of the metric, equation fIB.Qp shows 
that we are dealing with a generalised "cone" whose base is given by two spheres, §3 x §2 
|47] . However, the defining equation of the conifold is more readily written in terms of 
four complex coordinates Wi, i = 1 . . . 4 which obey [97] 
4 

i=l 

From equation (IB. 1211 . one sees that the geometry becomes singular at = for all 
i = 1 . . .4. In this case, §2 and S3 are shrinking to a point at the tip of the conifold 
where r = 0. Using the six- dimensional metric ( IB. 91) within the ten-dimensional ansatz 
f lB.7p . the components of the Einstein equations can be combined to yield the differential 
equation 

^h] + ^-(^h]=0, (B.13) 



whose solutions read 

Mr) = C2 + ^. (B.14) 

Note that it is possible to choose C2 = and that the ten-dimensional space is Ricci 
flat. Equation flB.lSp can indeed be obtained by setting the Ricci scalar i? = 0. 

For the time being, we have considered the case of a simple (and singular) conifold, 
where a priori no background fluxes have to be present, i.e. it is consistent to choose 



Brane inflation and the WMAP data 



57 



B(2) = C(2) = 0. One can now turn these fluxes on; the resulting effect is to induce 
corrections in equation flB.14|) which are proportional to the ratio of the flux quanta 
Ad, /C in the NS-NS and the RR sector respectively, namely (d5(2)) and (dC(2)). This 
was the route taken in references [HI HQ]. In the presence of fluxes, the singular conifold 
is replaced by its deformed counterpart and this amounts to modify equation (1B.12I1 
into 

4 

^^f = z, (B.15) 

i=l 

where the warping number 2; is a non- vanishing complex number. The solution for h{r) 
in the deformed conifold case (again found from the Einstein equations or from i? = 0) 
is [51] 



h{r) = C3 + C4 



1 1 , A '^z' 

In 1 + 



(B.16) 



ISr^z"^ 162^4 \^ ^2 

The constants C3, C4 can again be determined from the chosen compactifled geometry. 

In "cone" terminology, equation flB.lSP implies that the §2 still shrinks to a point 
at the tip of the cone, but the S3 remains flnite. As mentioned above, this is achieved 
by turning on background fluxes for (d-B(2)) and (dC(2)), i.e. aligning them with two 
Poincare dual 3-cycles. Two such cycles which verify equation (]B.12|) are described by 
4 3 
A:Y,^] = z, 'B:xl-J2y' = z, (B.17) 

i=l i=l 

where the complex coordinates Wi have been decomposed into their real and imaginary 
parts, Xi and yi respectively. Alignment of the 3-forms leads to a Dirac quantisation of 
the fluxes along the 3-cycles 

^ I (dC„) =2.M, ^ I (dB,,) = -2.K . (B.18) 

These conditions impose restrictions on the possible ansatze for the flelds, and eventually 
on the constants C3, C4 in equation flB.161) . 

The difference between the simple and the deformed conifold will be most 
remarkable close to the tip of the cone at r = tq. Let us notice that tq ^ since 
by deforming the conifold, we have replaced the "sharp" tip by a flnite size sphere S3. 
However, at some distance from the tip, both metrics should essentially give the same 
description of the geometry and in particular, we can "glue" the warped metric flB.16P 
to its simple version ( ]B.14[) at a certain radius r^,. From this matching, the constants 
Ci and C2 become some function of C3, C4 and r^,. Far from the singularity, the metric 
(IB. 141) expressed in terms of those constants provides an accurate enough description 
of the warped conifold geometry. Let us now link this setup with the brane inflation 
model we are interested in. We place two branes (D3 and anti-D3) into the background, 
and r will denote the coordinate distance between them. Now, one of these branes (the 
anti-D3 brane) is placed at the tip rg while ensuring r to be large. In that case, the 
metric in the vicinity of the D3 brane located at ri S> vq can still be accurately described 



Brane inflation and the WMAP data 



58 



by equation (lB.14p . It is for this reason that, while the solution (1B.16P for the warped 
metric is known, for all "practical" applications the simple solution flB.141) suffices. 

While we do not write down the exact solutions for the other fields in this theory 
exactly, let us state that the solution for the 4-form C(4), in the string frame, is found 
to be [221 iHl SHI Em El [97] 

[^(4)]oi23W=^s"^^- (B.19) 

This form turns out to be crucial in the derivation of the effective inflaton potential as 
it is demonstrated below. 



Appendix B.2.2. Interpretation of the background fluxes. We conclude our investigation 
of the background with some comments on an intuitive interpretation of background 
fluxes. If we consider a simple (singular) conifold geometry with non-zero fluxes, it can 
be viewed as a stack oiM = M.1C branes with RR charge placed at the r = singularity. 
The total RR charge N is therefore the product of the fluxes. If one uses the deformed 
conifold instead, this singularity is replaced by a three-sphere §3 with radius ro, and the 
charge M of the branes is literally smeared out across the three-sphere. In the following 
sections, we flrst add one more D3 brane (so that there are {M + 1) RR charged branes 
in total) to this background. This additional brane can be treated as a perturbation of 
the metric background ( IB. 141) . We then probe the perturbed background using a test 
anti-D3 brane and calculate the force it experiences. 



Appendix B.3. Perturbed geometry 

In order to describe a Dp brane inserted into this super-gravity background, one has to 
consider two additional contributions to the action (IB.ip [57j. In the string frame, they 
are the DBI action of a Dp brane, 

Sdp = -/Up / d^+i^e-* [- det (G„^ + B^^ + 27^a' F^^)f'^ , (B.20) 
plus its Chern-Simons coupling 

Sc-S,branc = ±Q'p / ^ (p+l) ! (B-21) 
J Mp 

since Dn branes in type IIB Superstring Theory are localised sources of C(p+i) fleld 
strengtlo. The RR charge is given by = jjip. For stabilised dilaton = e^*\ the 
effective brane tension Tp is related to the fundamental brane tension /ip by 

Tp = ^, (B.22) 

and can be calculated from [57] 

Tl = — (47rV)'~^ (B.23) 

* More generally, a Dp brane in type IIB theory can source C(;) fields for all even I <p-\-l [57] . 
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The gravitational coupling constant kiq was given in equation The upper sign 

in equation flB.21|) refers to the case of a Dp brane, the lower to an anti-Dp brane. 
The quantities ^* are the coordinates on the brane, Gaf3 is the induced four-dimensional 
metric and B^f^ the induced anti-symmetric tensor (Kalb-Ramond field) on the brane. 
These two are given by 

Gap = 9mn-^-Q^ , (B.24) 
Bap = Bmn-g^g^ , (B.25) 

gmn and Bmn being ten-dimensional. It is moreover possible (notably for parallel branes) 
to choose a configuration in which the gauge field on the brane Fap vanishes. 

In our case, the brane has three spatial dimensions, and we may choose the ^° to 
be parallel to the x^), i = 1, 2, 3 coordinate axes. If the coordinates y^, A = 4 . . .9 
of the brane are assumed to be static and independent of the location on the brane, it 
can be shown that there is no 5(2) "pull-back" onto the brane. Then equation (]B.24p 
implies 

Goo = goo, Gii = gu, for i = 1,2, 3. (B.26) 

The determinant of this induced metric, given that the ten-dimensional metric is (lB.7p . 
reads 

detG„;3 = G' = --^. (B.27) 

Notice that we would pick up an additional a^{t) factor if we were to use a FLRW 
universe on the brane instead of rj^i,. 

D3 branes are sources of C(4) fields, and thus of the corresponding F(5) field strength, 
as can be seen from the Chern-Simons coupling (IB. 211) 

dx™ dx" dx'^ dx' 



the summation over i running over all the D3 branes present ITU] . In this equation, 
the integration is carried out over the world volume M^^ of each brane, and /ig denotes 
their respective RR charge. It then follows from equations flB.20p and flB.2ip that the 
presence of one D3 brane at ri adds the following terms to the action 

-Ts / d'^V^ + fis I Ci^,), (B.29) 

J Ma J Ma 

where G is given by equation flB.271) . The D3 brane will manifest itself as a small 
perturbation of the background, i.e. the function h{r) calculated in Appendix B.l is 
replaced by 

h{r,ro) =h{r)+5h{ro). (B.30) 

Let us consider the effect of replacing h{r) by h{r,ri) in the Einstein equations. Since 
h{r) already fulfils equation (]B.13p . the resulting equation for 6h{ri) is 

^2 .\ 5 / 9 . \ _ _ , 



CO / j^j^Wj^(i) 



JDS 
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which has the solution 



Shin) =c, + ^. 



Arf 



(B.32) 



The new constants Ci and C2 can be related to the former Ci and C2 in equation (IB. 140 
by noting that the difference between the previous scenario and the present case comes 
from adding one D3 brane. One may therefore intuitively guess that [22] 



for z = 1, 2 , 



(B.33) 



The perturbed background metric h{r, ri) is obtained from equations flB.141) and flB.32p 
and, without loss of generality, one can set C2 = C2 = to get 



4^4 



(B.34) 



Note that as a consequence of the presence of the D3 brane (which is a localised RR 
charge source), the C(4) field is also perturbed, as can easily be seen from fIB.lQp . 



Appendix B.4. Test anti-brane in the perturbed geometry 

Let us denote by |/q = tq the radial position at which we insert the anti-D3 brane, while 
for all other extra coordinates we take Uq = y^, A = 5 ... 9. The anti-brane has the 
same tension but opposite charge as the brane and its contribution to the action in the 
string frame reads 

§D-3 = -^3 / d^e - /is / C(4) , (B.35) 

G' being the induced four-dimensional metric on the anti-brane. We assume that the 
anti-brane is "light" in the sense that it does not affect the already perturbed geometry. 
While we keep the ?/q static [i varying from to 5), we assume that the radial anti- 
brane coordinate y^ = tq does depend on the time x° (but not on the x*). Deriving G' 
explicitly leads to the expression 

2" 



det Gi 



al3 



G' 



1 



/i2(ro,ri) 



1 - /i(ro,ri) ( ^ 



(B.36) 



Inserted into the perturbed geometry of the previous subsection, i.e. where the warp 
factor is /i(r, ri) and the 4-form C(4) is also perturbed, the mobile anti-brane contributes 
to the action the terms 

1/2- 



-•ds 



/i(ro,ri) 



1 + 



(B.37) 



where we have used equation (]B.19|) . Expanding this expression for small velocities 
yields 



d^e 



2 



drg 
dt 



2T3 
h{r) 



(B.38) 
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By applying the canonic renormalisation = ^/^3r, the first term in (1B.38P becomes 
a standard kinetic term for a scalar field 0. The second term in equation flB.38|) can 
then be understood as the field potential and has to be evaluated at the position of the 
anti-brane. By means of equation (lB.34p . one gets 

^(-0, n) = = . (B.39) 

In the limit where the branes are far apart, namely ri -C ro, one has 

It can moreover be shown that Ci is related to M and v through 

Ci = 167r^(?sa" = 4r^v, (B.41) 

where r^^ is the "edge" of the KS throat (see section [23!) . 

In this Appendix, we have recalled the derivation of the interaction potential 
between the brane and the anti-brane in a ten-dimensional background warped by fluxes. 
The resulting potential would be the same if one chose to insert the anti-brane first into 
the background, and then study the motion of the brane. The only difference is that 
a D3 brane, inserted into the warped background, does not feel any force because of 
the compensation between the Chern-Simons coupling and the gravitation due to its 
tension [note the opposite signs of the two terms in equation flB.29P ]. This is not the 
case for an anti-brane inserted first, for which the two forces add up due the same signs 
in (IB. 351) . Seeking to minimise its energy, the anti-brane therefore sinks to the bottom 
of the KS throat. Let us recall that the above treatment remains valid only if ri ^ rg 
and breaks down as soon as the two branes become too close (see section [231). 
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